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MODELE^G  OF  OCEAN  SURFACE  SPECTRUM  AND  SURFACE 
SCATTERING  —  APPLICATIONS  TO  REVERBERATION 
PREDICTION  IN  THE  AAUS  PROGRAM 


1 .  Introduction 

There  are  two  approaches  to  modeling  long-range  reverberation  in  under¬ 
water  acoustics:  one  is  to  describe  reverberation  by  a  multi-element  model, 
of  which  each  element  separately  treats  surface  scattering,  bottom  scattering, 
propagation  within  the  medium,  and  complex  finite  target  effects;  and  the 
other  is  to  construct  a  reverberation  model  essentially  as  a  propagation  model 
in  which  surface,  bottom  scattering,  and  complex  target  scattering  are  taken 
into  account  via  application  of  appropriate  boundary  conditions.  The  modular 
approach  offers  a  practical  flexibility  that  allows  one  to  take  full  advantage 
of  a  broad  range  of  research  advances  attained  in  diversified  research  areas 
that  are  normally  investigated  independently  of  each  other.  For  this  reason 
the  NRL  reverberation  model  (RASP:  Range-dependent  Active  Systems  Prediction 
Model)  is  a  multi-element  model.  (For  a  documentation  of  RASP,  please  refer 
to  NRL  Report  No.  8721,  "NRL  Reverberation  Model:  A  Computer  Program  for 
the  Prediction  and  Analysis  of  Medium  to  Long-Range  Boundary  Reverberation," 
by  E.  R.  Franchi,  J.  M.  Griffin,  and  B.  J.  King,  1984.) 

There  have  been  several  methods  of  formulating  the  surface  scattering 
problem.  Some  were  developed  further  than  others,  but  each  has  its  own  limita¬ 
tions  and  strengths.  In  this  report,  these  different  methods  will  be  identified 
either  by  their  solution  methodology  or  by  their  fundamental  assumptions  if 
the  same  analytical  framework  is  to  be  taken.  The  two  major  methods  that  have 
been  subject  to  more  extensive  investigation  are  the  perturbation  and  tangent 
plane  methods.  In  addition  to  these  two,  other  methods  that  are  built  on  the 
concepts  of  composite-roughness  and  empirical  modeling  will  be  reviewed  and 
assessed  in  this  report.  The  current  effort  is  thus  focused  upon  those  surface 
scattering  models  that  deal  with  a  randomly  rough  surface  and  that  are  relevant 
to  the  modular  approach.  Therefore,  excluded  from  discussion  in  this  study 
are  research  works  [J1-J6]  treating  surface  scattering  as  part  of  a  propagation 
model  (e.g.  coupled-mode  analysis,  transport  theory)  and  those  developed  for  a 
periodically  rough  surface. 

While  RASP  is  a  complete  executive  routine  and  model  suite  capable  of 
handling  range  dependence,  complex  surface  and  bottom  geometries,  and  target 
characteristics,  its  operation  has  been  slowed  by  separate  calculations 
required  for  surface  scattering.  To  better  organize  and  thereby  facilitate 
the  application  to  extensive  data  bases  and  system  performance  calculations 
the  direct  incorporation  of  several  surface  scattering  submodels  has  been  under¬ 
taken.  The  program  has  entailed  a  critical  review  and  assessment  of  existing 
models  for  ocean  surface  scattering,  selection  of  several  surface  scattering 
models  that  are  deemed  most  suitable  and  incorporation  of  these  as  optional 
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submodels  addressable  by  the  RASP  executive.  Selection  of  the  scattering 
models  targeted  for  detailed  examination  and  earliest  incorporation  has  been 
guided  by  the  presently  prevalent  interests  in  low  frequencies,  low  grazing 
angles,  and  long  ranges  of  propagation. 

Despite  numerous  differences  among  these  surface  scattering  models,  predic¬ 
tions  made  by  them  are  all  critically  dependent  on  an  adequate  understanding 
of  the  ocean  surface  waves.  In  fact,  ocean  surface  wave  spectra  as  the  most 
important  environmental  parameter  strongly  influence  the  computational  results 
generated  by  these  models.  Comparison  with  experimental  measurements  and 
among  themselves  provide  resolution  only  when  ocean  surface  wave  spectra  are 
reasonably  understood  and  specified.  Consequently,  a  comprehensive  treatment 
of  the  ocean  surface  wave  spectrum  is  also  an  essential  requirement  in  making 
the  best  use  of  existing  scattering  theories  as  well  as  in  properly  including 
the  environmental  conditions.  Modeling  the  ocean  surface  wave  spectrum  has 
therefore  been  considered  an  integral  part  of  the  current  program. 

Currently,  a  comprehensive  review  of  surface  scattering  models  has  been 
completed.  These  scattering  theories  have  been  screened  and  assessed  with 
regard  to  their  relevance  to  current  interests.  The  scattering  models  that 
appear  most  practically  implementable  received  detailed  examination  in  this 
report.  The  existing  one-  and  two-dimensional  models  of  the  ocean  surface 
waves  have  also  been  reviewed  and  model  subroutines  for  these  surface  wave 
spectra  have  been  established  as  necessary  secondary  support  for  the  scattering 
models.  Preliminary  application  of  the  developed  one-  and  two-dimensional 
models  for  the  ocean  surface  wave  spectrum  has  been  Implemented  for  several 
chosen  scattering  models. 

The  general  purpose  of  this  report  Is  to  document  and  summarize  the 
research  progress  made  to  date  on  implementing  surface  scattering  acoustic 
models  and  their  component  ocean  surface  models,  as  components  of  RASP.  It 
will  also  chart  out  the  future  research  plan  for  this  effort.  Section  2  will 
provide  an  overview  of  major  methodologies  that  have  been  applied  to  investigate 
the  problem  of  surface  scattering.  A  detailed  examination  of  the  more  estab¬ 
lished  scattering  models  and  ocean  surface  spectrum  models  adapted  for  use  in 
the  selected  scattering  models  will  be  presented  as  appendices  to  the  main  text. 
Section  3  then  will  assess  the  relevance  of  each  major  scattering  model.  In 
Section  A,  computational  examples  of  several  scattering  models  combined  with 
the  chosen  ocean  surface  descriptions  will  be  Illustrated.  The  final  section. 
Section  5,  will  summarize  the  key  conclusions  and  planned  future  efforts. 

Also  provided  is  a  selected  bibliography,  in  which  key  references  of  the  same 
research  focus  are  grouped  together  at  the  authors'  discretion. 
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1 .  Alternative  Methodologies  in  Surface  Scattering  Modeling 

The  physical  picture  of  scattering  and  reflection  of  acoustic  waves  from 
the  ocean  surface  is  a  stochastically  random  process,  both  in  space  and  time. 
Fundamentally,  analyses  of  this  problem  can  be  formulated  from  either  the 
differential-equation  or  integral-equation  statement  of  the  problem.  From  the 
integral-equation  formulation,  one  can  single  out  the  effects  on  the  sound 
field  due  to  the  ocean  surface;  while  the  solution  to  the  differential- 
equation  statement  of  the  problem,  in  principle,  is  more  complete  but  it  is 
less  separable  since  scattering  and/or  reflection  is  included  as  part  of  a 
boundary-dependent  propagation  solution.  As  implied  from  the  stated  purposes 
of  this  study,  subsequent  discussions  will  be  primarily  concerned  with  the 
integral-equation  framework. 

In  evaluating  different  methodologies  of  solving  the  surface  scattering 
problem,  a  general  consideration  involves  identification  of  the  type  of 
boundary  description  for  the  ocean  surface.  Exact  analytical  solutions  are 
possible  when  one  treats  the  ocean  surface  as  a  periodically  rough  surface. 

The  most  general  statement  of  this  problem  specifies  the  ocean  surface  as  a 
completely  stochastic  ensemble  of  propagating  sinusoidal  waves  spanning  the 
entire  frequency  spectrum.  In  some  cases,  the  problem  is  simplified  somewhat 
by  eliminating  the  time  dependence  and  assuming  the  ocean  surface  as  a  ran¬ 
domly  rough  frozen  surface.  Excluding  the  case  of  periodically  rough  surfaces, 
this  study  will  consider  the  randomly  rough  surface  and  will  specify  whether 
It  Is  a  static  or  dynamic  one  for  the  methodology  under  consideration. 

In  addition  to  the  boundary  description,  other  general  considerations  that 
are  needed  In  the  process  of  assessing  alternative  methodologies,  are  associated 
with  the  chosen  geometry  of  source  and  receiver  as  well  as  the  fundamental 
characteristics  of  the  water  medium.  Scattering  solutions  reported  In  the 
literature  are  normally  tied  to  a  specific  source/receiver  geometry.  Since 
the  Intention  Is  to  Incorporate  a  surface  scattering  submodel  Into  a  general 
procedure,  this  study  will  focus  on  the  scattering  cross-section  solution, 
which  describes  the  scattering  Intensity  In  each  scattering  direction  per  unit 
area  per  unit  solid  angle.  For  the  same  reason  as  well  as  the  above-mentioned 
emphasis  on  the  Integral-equation  formulation,  effects  due  to  propagation-related 
medium  features  will  be  considered  In  the  total  calculation,  evaluated  In  the 
framework  of  RASP. 

However,  the  most  Important  consideration  that  helps  to  differentiate 
various  methodologies  rests  with  the  fundamental  approximation  assumptions 
employed  by  each  methodology.  Implicitly  Implied  from  these  assumptions  will  be 
the  validity  ranges  of  the  solutions  derived  from  the  corresponding  methodologies 
It  Is  therefore  most  natural  to  classify  existing  numerous  solutions  to  the 
surface  scattering  problem  principally  on  the  basis  of  their  fundamental  approxi¬ 
mation  assumptions  In  combination  with  their  solution  approach.  Discussions  on 
various  approximations  are  considerably  facilitated  through  the  use  of  several 
characteristic  lengths  associated  with  the  problem.  They  can  be  identified  as: 

a,  the  rms  surface  wavehelght; 

%  ,  the  acoustic  wavelength;  and 
A  ,  the  ocean  surface  wavelength. 
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Even  though  one  may  assume  Chat  these  approximations  can  all  be  trans¬ 
lated  into  some  well-defined  mathematical  statements  involving  only  these 
three  characteristic  lengths,  the  task  of  examining  various  approximation 
schemes  is  yet  further  complicated  by  the  fact  that  the  same  approximations 
are  usually  stated  in  different  forms  or  implemented  in  different  fashions. 
Therefore,  at  the  expense  of  neglecting  some  of  Che  more  subtle  aspects,  the 
following  discussions  on  alternative  methodologies  will  be  carried  out  in 
light  of  the  two  most  often-used  fundamental  approximation  assumptions.  These 
are  the  small-waveheight  and  Kirchhoff  approximations.  They  involve  specific 
relationships  between  the  above-mentioned  characteristic  lengths. 

It  should  be  noted  at  this  point  that  there  exists  another  fundamental 
approximation  assumption  that  is  widely  made  with  regard  to  the  time  variation 
of  the  ocean  surface.  According  to  this  assumption,  the  time  variation  of  the 
ocean  surface  is  relatively  slow  in  comparison  with  the  velocity  of  wave  prop¬ 
agation.  This  quasistatic  picture  allows  for  a  mathematically  more  manageable 
but  reasonable  description  of  the  time-dependent  scattered  field.  In  dealing 
with  the  dynamic  ocean  surface,  one  would  expect  some  variations  in  the  spec¬ 
trum  of  the  scattered  field  relative  to  that  of  the  incident  field.  This 
assumption  then  leads  to  an  a  priori  observation  that  the  scattered  acoustic 
field  will  be  a  narrowband  signal  centered  about  the  source  frequency.  As  a 
result,  this  quasistatic  approximation  of  the  dynamic  ocean  surface  is  nor¬ 
mally  referred  to  as  the  narrow-band  approximation. 

The  small-waveheight  approximation  is  most  often  applied  to  the  case 
where  the  slope  and  amplitude  of  the  surface  roughness  are  small.  By  this  it 
is  meant  that  the  rms  waveheight,  a,  is  small  compared  with  the  acoustic  wave¬ 
length  A.  and  the  minimum  ocean  surface  wavelength  Aniin  (  -A  niin  is  associ¬ 
ated  with  the  maximum  cut-off  frequency  of  the  ocean  surface  spectrum  being 
considered).  Mathematically,  one  writes  this  approximation  as 

a  «  \  and  a  «  /V  .  (2.1) 

The  Kirchhoff  approximation  is  usually  invoked  when  the  scattered  field 
at  an  observation  point  is  written  in  the  form  of  a  Helmholtz  integral  using 
the  half-space  Green's  function.  This  necessitates  the  knowledge  of  the 
first  derivative  of  the  scattered  field  at  the  Interface,  a  quantity  that  is 
not  readily  available  in  the  case  of  a  randomly  rough  surface.  The  Kirchhoff 
approximation  assumption  then  states  that  the  required  derivative  is  equal  to 
the  first  derivative  of  the  incident  wave,  which  is  a  known  quantity.  The 
equivalence  of  this  assumption  is  the  statement  that  the  ocean  surface  is 
"locally  flat",  or  in  terms  of  the  chosen  characteristic  lengths, 

X  Amin  .  (2.2) 

In  examining  the  connection  between  the  small-waveheight  and  Kirchhoff 
approximations,  Lablanca  and  Harper  ICIJ  pointed  out  that  in  several  existing 
solutions,  the  Kirchhoff  approximafon  is  applied  in  conjunction  with  the 
small-waveheight  one.  In  such  cases,  the  validity  domain  is  specified  by 

a  «  X«  Amin  ;  (2.3) 
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and,  in  fact,  is  contained  within  an  overlapping  region,  that  represents  the 
limit  a/\ -*  0  for  the  Kirchtioff  approximation  and  A.  /  A-"*  0  for  the  small- 
waveheight  approximation.  The  domains  of  validity  of  these  two  approximation 
schemes  together  with  their  overlapping  region  are  illustrated  by  a  Venn  dia¬ 
gram  in  Fig.  1.  In  this  report,  the  major  methodologies  are  examined  in 
Appendices  A  through  D  with  regard  to  their  assumptions,  formulation,  and  fun¬ 
damental  results.  In  reference  to  Fig.  1 ,  one  can  readily  recognize  that  all 
perturbation  method  solutions  are  valid  in  region  A.  Nevertheless,  most  of 
the  considered  Kirchhoff  solutions  are  restricted  to  the  overlapping  region  C 
instead  of  region  B  and  their  utility  is  therefore  relatively  limited.  On 
the  other  hand,  the  composite-roughness  solutions,  in  principle,  is  valid 
throughout  the  entire  acoustic  frequency  span  and,  with  a  judicial  separation 
of  large-scale  and  small-scale  portions  of  the  ocean  surface,  it  should  be 
practically  valid  for  any  frequency  range  of  interest.  For  the  acoustic  fre¬ 
quencies  that  range  from  very  low  to  medium  values,  the  solutions  derived  by 
the  perturbation  method  and  composite-roughness  approach  are  probably  more 
appropriate  for  use. 
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3.  Relevance  to  ML  Reverberatior  Model  RASP 


In  the  preceding  sections  and  attached  Appendices,  we  have  summarized 
and  assessed  a  number  of  theoretical  solutions  to  the  problem  of  sound  scat¬ 
tering  from  ocean  surface.  It  is  our  desire  to  direct  our  immediate  efforts 
to  satisfy  the  need  of  reverberation  prediction  for  the  AAUS  program.  Conse¬ 
quently,  our  selection  of  scattering  theories  and  what  to  do  with  them  will 
be  guided  by  the  following  requirements:  first,  the  selected  scattering  sub¬ 
models  must  be  practically  implementable  and  preferably  compatible  with  the 
structure  of  the  current  NRL  reverberation  model,  RASP;  second,  the  result¬ 
ing  computational  tools  should  be  able  to  simulate  as  closely  as  possible  the 
actual  complex  environmental  conditions;  and  finally,  order  of  preference 
should  go  in  tandem  with  planned  field  experiments. 

The  first  requirement  implies  that  our  effort  should  be  focused  on 
single-scatter  models  which  are  readily  Incorporated  as  addressable  components 
of  the  RASP  executive  to  describe  elemental  encounters  with  the  ocean  sur¬ 
face.  The  second  calls  for  an  adequate  description  of  the  ocean  surface,  of 
which  some  forms  of  solution  can  take  advantage  of  but  others  cannot.  This 
feature  is  required  to  further  enhance  forecasting  capability  by  including 
annual  and  short  term  projections  based  upon  actual  environmental  conditions, 
spatially  varying  sound  velocity  data,  and  detailed  bathymetric  data.  Finally, 
the  last  requirement  suggests  that  scattering  models  operational  in  the  low- 
frequency  and  low-grazing  angle  ranges  will  be  preferred.  In  addition,  this 
requirement  also  demands  that  we  should  first  concentrate  on  investigation  of 
variation  in  surface  reverberation  due  to  changes  in  environmental  parameters, 
such  as  wind  speed  and  wind  direction,  before  studying  aspects  of  surface 
reverberation  such  as  frequency  spreading/shifting  and  temporal/spatial  coher¬ 
ence,  some  of  which  are  under  Investigation  elsewhere. 

In  view  of  the  above  guidelines,  the  tangent  plane  method  solutions  exam¬ 
ined  in  Appendix  A  appear  not  to  address  the  frequency  range  of  interest  due 
to  their  underlining  Kirchhoff  assumption,  which  requires  the  acoustic  wave¬ 
length  be  much  smaller  than  the  surface  wavelength.  In  anticipation  of  this 
implication  of  the  Kirchhoff  approximation,  we  have  specifically  examined  the 
low-frequency  versions  of  these  tangent  plane  method  solutions,  which  were 
called  "long-wave"  by  Eckart  and  "slightly-rough"  by  Parkins.  Reduction  to 
this  special  case.  Indeed,  was  achieved  via  an  additional  approximation  that 
the  surface  wave  height  is  much  smaller  than  the  acoustic  wavelength,  not  to 
mention  the  introduction  of  the  small-slope  assumption  into  Eckart 's  solu¬ 
tion.  As  a  result,  the  use  of  these  additional  approximations  in  this  case 
ultimately  narrows  considerably  their  validity  domains.  A  more  rigorous  solu¬ 
tion  which  was  also  based  on  the  Kirchhoff  approximation  was  developed  by 
Tung,  who  did  not  use  the  small-slope  assumption  and  replaced  the  Fraunhoffer 
approximation  with  the  Fresnel  one.  His  solution  has  led  to  many  significant 
physical  insights,  especially  those  on  the  spectral  distribution  of  the  scat¬ 
tered  field.  However,  quantitatively,  the  form  of  his  solution  was  derived 
specifically  for  a  rather  simplistic  analytical  model  of  the  surface  wave 
spectrum  and  generally  cannot  accept  as  input  a  more  practical  surface  spec¬ 
trum  model  as  required  in  the  present  study.  Because  of  the  above-mentioned 
reasons,  none  of  the  tangent  plane  method  solutions  examined  in  this  report 
appears  suitable. 
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In  Appendix  B  we  have  closely  examined  three  approaches  to  solving  the 
surface  scattering  problem  using  the  perturbation  method.  The  fundamental 
assumption  unde’’ .  inning  this  method  is  the  small-waveheight  approximation, 
which  implie"  ^hat  the  rms  waveheight  is  much  smaller  than  both  the  acoustic 
and  surface  wavelengths.  Their  domains  of  validity  therefore  coincide  with 
the  low-1 requency  range  of  our  interest.  A  common  feature  of  these  solutions 
is  that  they  all  consist  of  two  components,  one  describing  the  specular  reflection 
and  the  other  depicting  scattering  in  all  other  directions.  Most  authors  chose 
to  examine  in  detail  the  scattered  field  except  Kuperman  who  focused  mainly  on 
the  specularly  reflected  component.  The  choice  of  which  solution  to  be  used 
in  this  class  of  solutions  is  greatly  facilitated  due  to  the  fact  that  the 
solutions  examined,  although  derived  within  different  theoretical  formulations, 
agree  with  each  other  except  for  a  few  nonessential  details.  Furthermore,  we 
will  be  primarily  concerned  with  the  far-field  solutions,  which  are  what  we 
need  for  combining  with  a  ray-theoretic  framework.  In  view  of  the  above  obser¬ 
vations,  it  is  clear  that  a  well-established  solution  based  on  the  first-order 
perturbation  method  exists  and  would  satisfy  our  Interest  in  the  range  of  low 
frequencies.  Another  attractive  aspect  of  this  solution  is  that  it  provides  an 
adequate  description  of  surface  scattering  up  to  the  first  order  in  three  spatial 
dimensions  as  well  as  frequency  domain.  As  a  consequence,  the  perturbation 
method  solution  as  represented  by  Eqs.  (B.30)  and  (B.53)  is  recommended  for 
use  in  RASP. 

At  freqencies  greater  than  a  few  hundred  hertz  the  ratio  of  rms  surface 
height  to  acoustic  wavelength  can  become  significant  due  to  the  large  amplitude 
of  the  very  long  surface  waves.  The  Kirchoff  solution  is  adequate  for  the  des¬ 
cription  of  scattering  from  large  gently  undulating  waves  but  is  not  adequate 
for  the  treatment  of  the  diffractive  scattering  induced  by  the  small  amplitude 
short  wavelength  component  of  the  surface  roughness.  The  composite  roughness 
model  (Appendix  C)  treats  the  effect  of  the  small  scale  roughness  as  a  perturba¬ 
tion  of  the  Kirchoff  scattering  from  the  large  scale  roughness.  This  extension 
of  perturbation  theory  greatly  exends  the  domain  of  validity  of  the  perturbation 
approach.  The  theory's  most  serious  limitations  occur  at  very  low  grazing 
angles  where  the  Kirchoff  component  of  the  approximation  requires  corrections 
for  shadowing.  The  composite  roughness  submodel  is  summarized  by  Eqs.  (C.6) 
and  (C.28).  This  submodel  should  be  used  with  caution  so  that  the  rms  slope 
of  the  large  scale  roughness  of  the  specified  ocean  surface  is  not  on  the 
order  of  the  relevant  grazing  angles. 

Another  group  of  surface  scattering  descriptions  reviewed  in  Appendix  D 
were  constructed  by  the  empirical  modeling  approach  on  the  basis  of  available 
experimental  data.  They  represent  empirical  fittings  of  fairly  comprehensive 
sets  of  scattering/reverberation  acoustic  and  optic  data.  In  the  case  of  the 
Chapman-Ha rr is  and  Martin  backscattering  models,  their  empirical  formulae  have 
been  used  extensively  and  have  proved  to  be  relatively  accurate  within  their 
Intended  ranges  of  frequencies,  grazing  angles,  and  wind  speeds.  Their  utility 
is  nevertheless  quite  limited  since  they  only  describe  the  incoherent  scattering 
strength  for  the  backscattering  geometry.  Application  of  the  empirical  models 
thus  cannot  be  extended  to  general  scattering  geometries  and  spectral  analysis. 
However,  because  of  the  fact  that  they  were  worked  out  from  field  data,  they 
are  attractive  in  many  aspects  for  comparison  with  other  surface  scattering 
models.  RASP  already  includes  the  Chapman-Harris  model  and  the  Martin  model 
will  be  Incorporated  later. 
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So  far  our  selection  of  surface  scattering  submodels  has  included  solu¬ 
tions  by  the  first-order  perturbation  method  and  the  composite-roughness 
approach  and  empirical  formulas  proposed  by  Chapman-Harris  and  Martin.  There 
exist  other  physical  mechanisms  that  would  modify  the  scattered  field  but  were 
not  accounted  for  by  these  solutions.  Among  them,  the  two  most  important  ones 
are  bubble  scattering  at  low  grazing  angles  and  shadowing  effect  associated 
with  large-scale  roughness  and  low  grazing  angles.  Since  these  two  effects 
can  be  analyzed  separately  and  incorporated  additionally  to  any  surface  scat¬ 
tering  solutions,  we  intend  to  build  them  into  RASP  as  two  secondary  options 
which  can  be  combined  with  any  of  the  selected  surface  scattering  submodels. 

As  indicated  at  the  beginning  in  Section  1,  the  performance  quality  of 
any  surface  scattering  submodels  critically  depends  on  the  quality  of  the 
ocean  surface  description  in  use.  Specifically,  one  needs  to  have  as  inputs 
to  these  submodels  reliable  two-dimensional  spectra  of  the  surface  wave  ampli¬ 
tude  and  slope.  We  therefore  have  proceeded  to  review  major  modeling  efforts 
for  the  ocean  surface  waves  and  have  come  to  the  conclusion  that  the  best 
available  two-dimensional  surface  wave  spectrum  would  be  a  combination  of  the 
one-dimensional  spectrum  developed  by  Pierson  [K4]  and  the  azimuthal  distri¬ 
bution  proposed  by  Mitsuyasu  [K6].  This  combined  surface  wave  spectrum  is 
essential  for  use  in  the  surface  scattering  submodel.  In  Appendix  E,  this 
combined  model  is  reviewed  and  illustrated  with  sample  computed  spectra  of 
surface  waves. 

In  the  following  section,  the  recommended  solutions  of  scattering  cross- 
sections  will  be  computed  with  the  use  of  the  above  ocean  surface  model. 

Sample  computations  will  be  presented  to  illustrate  the  major  features  of 
these  solutions  and  their  variations  with  key  environmental  parameters. 
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4 .  Example  Computations 

On  the  basis  of  the  assessment  of  different  surface  scattering  solutions 
performed  in  the  preceding  sections  and  attached  Appendices,  we  have  reached 
the  recommendation  that  the  first-order  perturbation  method,  composite¬ 
roughness  approach  and  Chapman-Harris  empirical  model  be  used  for  computation 
of  the  surface  scattering  strength.  In  addition,  as  the  most  important  envi¬ 
ronmental  input,  the  two-dimensional  surface  spectrum  modeled  after  a  combine 
tion  of  the  six-range  Pierson's  one-dimensional  spectrua  and  the  Mitsuyasu 
azimuthal  dependence  function  will  be  employed  in  these  computations  whenever 
applicable.  The  objective  of  this  section  is  then  to  illlustrate  the  compu¬ 
tational  results  that  can  be  achieved  by  the  recommended  solutions. 


From  the  first-order  perturbation  method,  the  scattering  cross  section 
^ p[^  is  given  by  (see  Eq.  (E.30)  of  Appendix  E) 


(4.1) 


where  S2  is  Che  two-dimensional  surface  wave  spectrum.  Other  symbols  used  in 
Eq.  (4.1)  are  defined  as  follows: 

»  Incident  wave  vector, 
kg  =  Scattering  wave  vector, 

Wq  *  Carrier  frequency,  and 


tC  and  are  the  resonant  scattering  wave  vector  and  its  corresponding  Doppler 
shifted  frequency  to  be  specified  as 


K 


M  -  I 

*i.  > 


(4.2) 

(4.3) 

In  the  backscattering  case,  i.e.  kg  •*  -k^,  Eq.  (4.1)  is  reduced  to 


We  now  want  Co  examine  several  sample  computations  of  Che  scattering 
cross  section  with  the  use  of  Eq.  (4.4).  In  order  to  isolate  the  effect  of 
wind  speed  variation,  the  azimuthal  dependence  of  S2  is  removed  by  replacing 
it  with  an  isotropic  two-dimensional  surface  spectrum  S*  defined  by 


S  (K) 


1 

XT 


(4.5) 


where  Si  is  the  Pierson's  one-dimensional  spectrum  presented  in  Appendix  E. 
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Figures  2a  and  2b  present  the  variation  of  the  isotropic  backscattering 
cross  section  with  respect  to  representative  values  of  the  wind  speed  at  two 
acoustic  frequencies,  200  Hz  and  50  Hz,  respectively.  In  the  200  Hz  case, 
the  variation  of  (T with  wind  speed  appears  minimal  (less  than  5  dB)  for  all 
grazing  angles,  except  for  those  angles  higher  than  60°.  This  is  mainly 
influenced  by  the  fact  that  at  200  Hz  the  resonant  scattering  frequency. 


varies  from  0.65  Hz  to  0.25  Hz  as  the  grazing  angle  increases  from  0°  to  70°. 
In  this  frequency  range  of  surface  waves,  it  is  seen  from  Fig.  E2  that  the 
spectral  level  does  not  vary  considerably  with  the  wind  speed.  At  the  lower 
acoustic  frequency,  50  Hz,  Fig.  2b  shows  that  that  variation  of  the  isotropic 
backscattering  cross  section  with  respect  to  wind  speeds  is  more  significant, 
especially  between  3.5  m/s  and  7.4  m/s.  This  characteristic  again  can  be 
related  to  the  surface  spectrum,  where  the  spectral  level  is  seen  to  vary  con¬ 
siderably  in  the  frequency  window  0.12  -  0.32  Hz,  representing  the  variation 
of  the  resonant  scattering  frequencies  in  the  same  grazing  angle  range. 


With  regard  to  the  composite-roughness  approach,  the  general  functional 
behavior  of  the  scattering  cross  section  is  similar  to  that  of  the  perturba-  . 
tion  method,  except  for  the  tilting  effect  that  modifies  the  effective  grazing 
angle.  Its  scattering  cross  section  can  be  written  as  (see  Eqs.  (C.6)  and 
CC.28)) 

■  ^(-X)SCto-u3,+x^)l  , 

in  which 

s 

e 

where  and  Og  are  the  grazing  and  scattering  angles,  €  the  rms  slope  of 
surface  waves,  and  other  symbols  are  defined  as  before.  The  above-mentioned 
tilting  effect  is  embodied  in  the  quantity  of  Eq.  (4.7). 


=  -  ^  ^  f !  V  ^ C  ( V  ^  ^  3 


=  ©  1-  e, 


S  ) 


= 


(4.8) 


s  ; 


In  the  case  of  backscattering,  Eqs.  (4.7)  and  (4.8)  are  simplified  to 


in  which 


Again,  we  compute  and  present  in  Fig.  3  the  isotropic  backscattering 
cross  section  as  derived  by  the  composite-roughness  approach  for  different 
wind  speeds  at  the  acoustic  frequency  of  200  Hz.  Variation  with  the  wind 
speed  is  seen  not  to  be  considerable.  In  this  particular  set  of  calculations, 
we  also  include  a  computation  that  uses  the  familiar  Phillips  spectrum  as 
input  instead  of  our  adapted  Pierson's  spectrum.  The  resulting  difference  is 
obvious  as  indicated  in  Fig.  2,  wtiich  demonstrates  the  criticality  of  which 
surface  spectrum  model  to  be  used  as  input  in  theoretical  prediction  of  the 
surface  scattering  cross  section. 

The  above  two  surface  scattering  theories  together  with  the  Chapman- 
Harris  empirical  model  (see  Appendix  D)  are  now  compared  using  the  same  acous¬ 
tic  and  environmental  parameters.  Figure  4  presents  the  isotropic  backscat¬ 
tering  cross  section  computed  by  these  thre>  models  for  acoustic  frequency  of 
200  Hz  and  wind  speed  of  3.5  m/s.  At  high  grazing  angles,  e.g.  larger  than 
50°,  their  predictions  are  very  similar.  However,  as  the  grazing  angle 
becomes  small,  their  predictions  begin  to  deviate  from  each  other.  Chapman- 
Harris  model's  prediction  is  about  at  least  10  dB  consistently  lower  than 
those  of  the  other  two  solutions.  Among  the  perturbation  method  and  composite¬ 
roughness  approach,  the  scattering  strength  predicted  by  the  latter  approaches 
a  low  but  finite  level  at  about  -70  dB  for  grazing  angles  of  a  few  degrees, 
while  the  former  practically  predicts  no  backscattered  energy.  In  general, 
the  difference  between  the  perturbation  method  and  composite-roughness 
approach  is  only  significant  for  grazing  angles  less  than  20°. 

In  the  above  computations,  we  have  removed  the  azimuthal  dependence  of 
the  surface  scattering  cross  section  and,  as  a  consequence,  effectively  inco¬ 
herently  summed  the  two  constituent  Doppler-shifted  components.  Utilization 
of  the  proper  two-dimensional  surface  spectrum  instead  of  Eq.  (4.5)  would 
enable  us  to  resolve  these  two  components  and  examine  the  effect  of  the  wind 
direction  upon  their  respective  strengths.  In  Figs.  5a  through  5d,  the  two- 
dimensional  surface  spectrum  developed  in  Appendix  E  is  employed  to  evaluate 
the  backscattering  strength  of  both  the  up-shifted  and  down-shifted  frequency 
components  for  different  wind  directions  but  the  same  wind  speed.  As  the  wind 
direction  is  changed  from  almost  cross-wind  (Fig.  5a)  to  up-wind  (Fig.  5d), 
the  down-shifted  component  slowly  gains  strength  while  the  up-shifted  compo¬ 
nent  quickly  loses  its  intensity.  For  the  cross-wind  condition,  the  two  com¬ 
ponents  are  expected  to  be  equal  in  strength  and  for  the  up-wind  condition 
only  the  down-shfited  component  exists. 

Indeed,  other  than  the  backscattering  case,  the  perturbation  method  as 
described  by  Eq.  (4.1)  can  predict  the  surface  scattering  cross  section  for 
all  scattering  directions,  both  in  and  out  of  the  plane  of  incidence.  In 
this  regard,  the  composite-roughness  solution  as  represented  by  Eq.  (4.7)  is 
able  to  describe  the  scattering  cross  section  for  all  directions  in  the  plane 
of  incidence.  These  capabilities  are  demonstrated  in  Figs.  6  and  7,  where 
the  incoherent  (up-shifted  plus  down-shifted)  scattering  strength  is  plotted 
versus  the  scattering  angle  in  the  plane  of  incidence  for  a  given  grazing 
angle.  The  computations  presented  in  Fig.  6  are  performed  for  wind  speed  of 
3.5  m/s,  acoustic  frequency  200  Hz,  grazing  angle  30°,  and  two  wind  directions 
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(up-wind  and  cross-wind),  employing  both  perturbation  and  composite-roughness 
solutions.  Similar  computations  are  presented  in  Fig.  7  using  the  same  para¬ 
metric  values  except  for  a  different  wind  speed,  7.4  m/s.  The  following 
observations  can  be  drawn  from  these  two  figures:  (i)  At  low  wind  speed, 

3.5  m/s,  the  scattering  strengths  vary  considerably  with  the  wind  direction; 
(ii)  At  moderate  wind  speed,  7.4  m/s,  and  at  acoustic  frequency  of  200  He, 
the  wind  direction  has  a  significant  effect  on  the  predicted  scattering 
strength  only  in  the  forward-scattered  directions;  and  (iii)  Difference 
between  the  perturbation  and  composite-roughness  solutions  is  noticeable  only 
for  scattering  angles  in  the  neighborhood  of  0®  and  180°. 


5 .  Conclusions  and  Future  Developments 

Assessment  of  various  theoretical  formulations  and  examination  of  their 
principal  results  in  this  report  has  led  to  an  appropriate  selection  of  alter¬ 
native  surface  scattering  submodels  and  associated  surface  wave  spectra  that 
would  be  most  compatible  with  AAUS  interests  as  well  as  the  current  NRL  rever¬ 
beration  prediction  model,  RASP.  The  requirements  of  low-  to  mid-range 
frequencies  and  low  grazing  angles  represent  the  driving  rationale  for  us  to 
adapt  the  following  approach  to  modeling  surface  scattering  in  our  current 
effort  in  reverberation  prediction: 

(i)  Two  theoretical  formulations,  the  first-order  perturbation 
method  and  the  composite-roughness  approach,  are  selected  to  be  incorporated 
into  RASP.  The  composite-roughness  solution  would  warrant  a  broader  cover¬ 
age  of  the  frequencies  of  interest.  The  principal  differences  between  these 
two  solutions  are  two-fold:  First,  near  0®  or  180®  scattering  angles  the 
composite-roughness  solution  approaches  a  finite  value  for  scattering  strength 
while  the  perturbation  method  leads  to  an  Infinitely  small  scattering  strength; 
and  second,  at  the  present  stage  the  composite-roughness  approach  can  describe 
the  scattering  cross  section  in  the  plane  of  incidence  whereas  the  perturba¬ 
tion  solution  is  capable  of  depicting  the  scattering  strength  in  all  three- 
dimensional  directions.  (One  reason  to  incorporate  the  perturbation  solution 
is  due  to  the  anticipation  of  the  eventual  consideration  of  three-dimensional 
effects  in  the  reverberation  problem.) 

(ii)  Two  backseat terlng  empirical  models,  Chapman-Harris  and 
Martin,  will  also  be  Introduced  into  RASP.  As  well-established  empirical 
formulae,  they  will  serve  as  useful  double-check  options  and  will  be  used  as 
bases  for  future  additions  of  new  field  measurements  into  an  empirical  fitting 
framework  as  needs  arise. 

(ill)  The  two-dimensional  surface  wave  spectrum  modeled  from  a  com¬ 
bination  of  the  six-range  Pierson  one-dimensional  spectrum  and  the  azimuthal 
distribution  proposed  by  Mitsuyasu  et  al.  appears  to  be  the  most  realistic 
and  comprehensive  empirical  description  based  on  available  information  at  the 
present  time.  In  addition,  the  choice  of  this  combined  spectrum  is  also  moti¬ 
vated  by  the  fact  that  it  adequately  meets  the  input  requirement  of  most 
existing  surface  scattering  formulations  and  its  fitting  algorithms  are  flex¬ 
ible  enough  to  accommodate  new  ocean  surface  measurements,  which  are  being 
acquired  in  the  present  experimental  plan  of  the  AAUS  program. 

At  the  completion  of  this  study,  our  effort  in  investigating  surface 
reverberation  has  by  no  means  reached  the  desired  understanding.  Future 
research  efforts  are  therefore  proposed  to  be  pursued  in  the  following 
directions: 

(1)  It  is  obvious  that  the  next  step  would  be  to  implement  the 
above  recommended  surface  scattering  options  in  RASP.  This  can  be  achieved 
with  some  minor  modifications  to  the  present  programming  code  of  RASP. 
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(ii)  Our  understanding  of  the  reverberation  intensity  would  be 
solidified  as  the  result  of  the  above  efforts.  However,  there  exists  a  strong 
need  to  also  characterize  the  spectral  content  of  the  surface  reverberation 
signal.  It  is  therefore  desirable  to  speed  up  investigation  of  the  problem 

of  frequency  shifting  and  spreading  for  surface  reverberation.  Due  to  the 
rather  immediate  need  of  this  understanding,  investigation  of  this  problem 
should  be  simultaneously  pursued  in  two  directions,  both  independently  of 
RASP  and  within  the  ray-theoretic  framework  of  RASP. 

(iii)  Finally,  in  a  more  long-term  time  frame,  the  stochastic 
fluctuations  of  the  ocean  surface  environment  which  would  obviously  degrade 
the  reverberation  signal  must  be  addressed  to  enhance  our  prediction  capabil¬ 
ities.  This  implies  that  we  will  have  to  investigate  the  amplitude  and  phase 
fluctuations  of  the  surface  scattering  field  and  consequently  can  estimate 
the  critical  concepts  of  temporal  and  spatial  coherence  of  the  surface  rever¬ 
beration  signal. 
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APPENDIX  A 


Tangent  Plane  Method 

In  general,  application  of  the  Kirchhoff  approximation  to  the  problem 
of  acoustic  scattering  from  the  ocean  surface  was  normally  done  within  the 
context  of  the  Helmholtz  integral  over  a  half-space  integration  contour. 

Eckart  lB7J  employed  this  formulation  to  investigate  scattering  from  a  frozen 
randomly  rough  surface.  Later,  Parkins  [85 J  extended  Eckart's  analysis  to 
include  the  time  variation  of  the  ocean  surface,  i.e.,  to  the  case  of  a 
dynamic  stochastic  surface.  Recently,  this  method  of  solution  was  further 
refined  by  Tung  lB4j  to  improve  its  accuracy  and  to  account  for  wave-slope 
statistics  and  a  more  general  source/receiver  geometry.  In  addition  to  these 
major  works,  several  other  authors  1C1-C3J  attempted  to  develop  ways  to  com¬ 
pensate  for  a  major  deficiency  Inherent  In  the  Kirchhoff  approximation, 
caused  by  the  so-called  shadowing  effect.  In  this  Appendix,  we  will  sum¬ 
marize  and  discuss  the  major  solutions  that  are  based  on  the  Kirchhoff  approx¬ 
imation  scheme.  Since  there  have  been  numerous  studies  using  this  approach, 
following  discussions  on  this  methodology  will  have  to  be  focused  on  the 
works  of  only  several  key  researchers.  Attempts  will  be  made  to  recast  these 
works  in  the  same  terminology  for  them  to  be  comparable. 

A1 .  Eckart 

The  problem  investigated  by  Eckart  lB7J  was  to  derive  the  scattering 
cross  section  for  an  acoustic  point  source.  This  was  formulated  from  the 
half-space  Helmholtz  integral,  which  describes  the  scattered  pressure  field 
as 


for  an  incident  field  of  the  form 


Utr. 


In  Eqs.  (A.l)  and  (A.2),  P  is  the  directional  source  amplitude;  ^  is 
the  ocean  surface  described  by  jbc  t};  tv  is  the  unit  vector  normal  to^; 

is  the  acoustic  wave  number;  and  and  ir  are  as  illustrated  in  Fig.  Al. 

At  this  point,  Eckart  applied  the  pressure-release  boundary  condition 
and  the  Kirchhoff  approximation  at  the  surface. 


=  0  , 

(A. 3) 

hi. 

3r, 

1 

(A. 4) 
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to  write  Eq.  (A.l)  as 


(A. 5) 


Equation  (A.l)  was  then  simplified  further  through  the  use  of  the  follow¬ 
ing  approximations; 


(i)  Fraunhoffer  approximation: 

r  r  R .  R  , 

(ii)  Small-slope  approximation: 

^  2. 

^  ^  ^  (A.7) 

“♦  J-kcL^  . 

The  Fraunhoffer  approximation  is  based  on  the  assumption  that  both  source 
and  receiver  are  in  the  far  field  and  the  source  beam  pattern  is  highly  direc¬ 
tional.  The  small-slope  approximation,  which  allows  the  integration  to  be 
carried  out  on  the  surface  the  projection  of  aC  onto  the  (x,y)  plane, 

implies  that  With  these  additional  approximations,  Eq.  (A. 5)  becomes. 


=  - 


i.U 


^7C 


(A. 8) 


A  A 

wher^ot ,  and  are  the  Cartesian  components  of(Rj-^^j.  If  we  define 
and  to  be  the  incid^t  and_^scattered  wave  vectors,  respectively,  then 
is  the  z-component  of  . 

4 


From  Eq.  (A. 8),  Eckart  proceeded  to  study  two  simplified  cases,  which 
were  called  "long-wave"  and  "short-wave"  in  his  terminology.  Since  the 
"short-wave"  solution  was  derived  for  a  specific  form  of  the  ocean  surface 
spectrum  and  thus  was  not  generally  useful,  we  will  discuss  only  the  "long¬ 
wave"  solution  here.  The  "long-wave"  case  was  prescribed  by  the  following 
approximation 


e 


I  -  MZ  , 


(A. 9) 


which  leads  to  the  condition 

2XZZ 
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or* 


«  X 


(A. 10) 
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Writing  the  intensity  of  the  scattered  f ield  from  Eq.  (A. 8)  in  terms  of 
the  spatial  covariance  function  of  the  surface  wave  height, 


(A. 11) 


Eckart  derived  the  scattering  cross  section  per  unit  solid  angle  per  unit 
area  as 


(T 

^c. 


,4TC 


(A. 12) 


where  ^()^,k^^is  the  spatial  two-dimensional  spectrum  of  the  ocean  surface 
waves,  which  is  defined  by 


s. 


(A. 13) 


for  a  homogeneous  ocean  surface 


A2.  Parkins 

The  problem  studied  by  Parkins  lB5J  differs  from  that  by  Eckart  in  two 
aspects:  An  incident  plane  wave  was  considered  instead  of  a  point  source 
and  time  variation  of  the  ocean  surface  was  included  in  Parkins'  treatment. 
Parkins'  solution  was  developed  from  a  half-space  Helmholtz  integral  of  a 
slightly  different  form,  which  was  based  on  an  impulse-type  Green  function 
Isee  J.  A.  Stratton,  "Electromagnetic  Theory,"  McGraw-Hill,  New  York,  1941, 
p.  427 J.  In  our  present  notation,  this  integral  can  be  written  as 


for  the  incident  field  given  by 

*  A0 

/•^  .  \  r=>  -  ^  uut: 

Ti  (•">)  -  ^  ^  ^ 


(A. 14) 


(A. 15) 


where  is  the  amplitude  of  the  incident  plane  wave  and  the  square  brackets 
denote  time  retardation.  Because  of  considering  ^n  Incident  plane  wave. 
Parkins  could  carry  out  the  differentiation  b/^fi  exactly  without  resorting 
to  the  small-slope  approximation  as  Eckart  did.  This  results  in 


(A. 16) 


in  which  the  Kirchhoff  approximation  has  been  applied  and  evaluation  at  the 
retarded  time  implemented. 
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Using  Eq.  (A. 16)  and  the  fact  that 


I  )  ) 


/a. 


(A. 17) 


Parkins  transformed  Eq.  (A. 14)  into 
p  _i«ot  ,, 


18) 


At  this  point,  Parkins  invoked  the  far-field  Fraunhoffer  approximation 
in  the  form 


r*  Cl  Rj  u.  XtM6)^  4uJl 

K  A  ■' 

CL  R5-  Rj.  r  =  Rj-  I'w  4U 


LlC^yyJ^ 


(A. 19) 


to  obtain 


“  '  **aic  R - f)  *  ^ 


Through  a  simple  Integration  by  parts,  Eq.  (A. 20)  becomes 

^  ,p  rr  .:K  r,  \  , 
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Now,  Parkins  made  the  assumption  that  "the  aperture  is  many  wavelengths 
across”  l85,  p.  1263J  to  obtain  the  following  key  result 


^  „  :Kt 


M5 


(A.22) 


From  Eq.  (A.22)  the  spectral  density  of  the  scattered  field  can  be 
formally  derived  from 


^  1P(^) 


-•0 


where  the  covariance  function  is  defined  by 


(A. 23) 


(A. 24) 
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The  averaging  operation  in  Eq.  (A. 24)  was  considerably  simplified  by 
assuming  that  the  joint  distribution  of  the  surface  wave  height  is  Gaussian. 
This  then  leads  to 


(A. 25) 


where 


^  £> 

"a- 


£  -  normalized  covariance  function  of  surface  waveheight, 

i 

rr  = 


The  spectral  density  of  the  scattered  field  was  then  solved  explicitly 
for  a  special  case,  which  was  called  "slightly  rough  surface"  by  Parkins. 
This  was  prescribed  by  the  approximation 


<C 


leading  to 


.  *■  *- 

~  k.T 


«  /  *  J.  _  \ 


(A. 26) 


Considering  the  insonified  area  as  a  square  of  (  it  )  dimension, 

Parkins  showed  from  Eqs.  (A. 23),  (A. 25),  and  (A. 26)  that  the  scattering  cross 

unit  area  for  the  "slightl^-^rough"  case  had 


section  per  unit  solid  aygle  per 
a  coherent  component  and  a  D 


Doppler-shifted  component  (p*  given  as 


r’u)  - 


Sc 


ATC 


(A. 27) 


(A. 28) 


where  is  the  two-dimensional  spectrum  of  surface  waves  defined 


by 
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A3.  Tung 


Tung  [B4J  generalized  Eckart's  solution  to  include  the  time  variation 
of  the  ocean  surface.  This  was  achieved  by  starting  also  from  Eq.  (A.l)  but 
the  scattering  pressure  fs  and  its  normal  direction  in  this  equation  was  evalu¬ 
ated  at  the  retarded  time.  A  major  improvement  in  terms  of  accuracy  of  Tung's 
solution  was  obtained  by  dropping  both  Fraunhoffer  and  small-slope  approxima¬ 
tions.  In  place  of  the  Fraunhoffer  approximation,  Tung  employed  the  more 
general  Fresnel  approximation.  Gain  in  accuracy  of  the  Fresnel  approximation 
is  attributed  to  the  fact  that  the  horizontal  components  are  kept  up  to  the 
second  order. 


This  can  be  better  illustrated  as  follows.  Referring  to  Fig.  A1 ,  one 
can  write  the  exact  expression  for  rg  as 


JjRf.r 


S-) 


(A. 30) 


from  which  the  Fraunhoffer  approximation  is  derived  in  the  form  (compare  with 
Eq.  (A. 19)) 


where  r  =  (x,y,Z).  More  generally,  if  one  retains  only  the  first-order  terms 
in  Z  but  keeps  the  x  and  y  terms  up  to  the  second  order,  one  obtains  the 
Fresnel  approximation.  From  Eq.  (A. 30),  this  approximation  is  described  as 


(A. 31) 


Tun^  also  avoided  the  use  of  the  small-slope  approximation  in  handling 
the  r\  differentiation  by  solving  the  Helmholtz  integral  for  a  specific 
surface  spectrum.  This  in  turn  considerably  limits  the  usefulness  of  Tung's 
solution,  especially  when  the  employed  form  of  the  surface  spectrum  is  fairly 
simplistic  and  completely  ignores  the  dispersive  character  of  ocean  waves. 
This  surface  spectrum  was  represented  by  the  covariance  function 


(A. 32) 


where  and  T,  are  the  correlation  lengths  with  respect  to  the  x,  y,  and 

t  dimensions.  Equation  (A. 32)  indeed  only  represents  a  surface  wave  propagat¬ 
ing  in  the  y  direction  with  a  velocity  C  and  wavenumber  K  . 
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With  regard  to  the  temporal  distribution  of  the  ocean  waves,  instead  of 
assuming  a  joint  Gaussian  distribution  like  Parkins,  Tung  carried  out  the 
averaging  operation  on  the  basis  of  the  assumption  that  the  surface  height  is 
"approximately"  Gaussian  and  is  uncorrelated  with  the  slope.  This  allowed 
Tung  to  establish  that  the  slope  statistics  was  one  of  several  reasons  that 
caused  asymmetric  sidebands  in  the  scattered  field. 

In  the  final  analysis,  Tung's  solution  is  more  rigorous  in  the  sense  that 
it  is  based  on  the  Kirchhoff  approximation  alone,  doing  away  with  the  far- 
field  Fraunhoffer  and  small-slope  approximations.  However,  its  application  is 
quite  limited  because  it  cannot  accept  a  general  empirical  or  experimentally- 
measured  surface  wave  spectrum. 


A4.  Shadowing  Function 

The  validity  of  various  solutions  to  the  surface  scattering  problem  based 
on  the  Kirchhoff  approximation  is  also  limited  by  not  incorporating  the  fact 
that  the  surface  is  only  partially  illuminated  due  to  the  roughness.  A  proper 
accounting  of  the  so-called  shadowing  effect  would  certainly  improve  the  rigor 
of  these  solutions.  The  fundamental  concepts  to  investigate  the  shadowing 
effect  was  first  developed  by  Beckmann  1G3J  and  further  elucidated  later  by 
Wagner  IGIJ.  According  to  Wagner,  the  probability  that  a  point  ,  with 
slope  Z{  of  the  rough  surface  is  illuminated,  can  be  solved  from  the  differ-  ' 
entlal  equation 


A[H(6j2„2»]  =  --l(X)H(6jZ,X;X)  ^  (A-33) 

in  which  the  spatial  parameters  are  illustrated  in  Fig.  A2,  W(6jZ,X,X)  is 
the  probability  that  the  surface  Z(X)  does  not  cross  the  ray  Z  =  Zi+XtanS^  in 
lO,Xj  and  q(X)d(X)  is  the  probability  that  Z(X)  crosses  Z  =  Zj+XtanSi  in 
lX,X+dXj  provided  Z(X)  does  not  cross  Z  =  Zj+XtanSi  in  10,XJ. 


The  solution  to  Eq.  (A. 33)  has  the  form 

H(e.|2„2,/)  =  axf  [.^  J  . 

0 


(A. 34) 


(A. 35) 


H(ej  %„  z,') »  i  ( ^  M  jlx]  ^ 
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where  h  is  the  unit  step  function  and  H(0i|Z]^,Z|)  is  the  probability  that 
the  point  Zj  with  slope  Z[  is  illuminated.  A  shadowing  function,  which  is 
defined  as  the  ratio  of  illuminated  to  total  areas  of  the  surface,  is  then 
obtained  by  averaging  over  all  possible  values  of  surface  height  and  slope, 

H  ( H  ( ft  jzX)  w  (z.)  I  (A.3 

— oA 

in  which  W(Zi)  and  W'(Zl)  are  the  probability  densities  of  surface  height 
and  slope,  respectively. 


At  this  point,  Wagner  employed  the  following  approximation  assumptions; 


(i)  q(X)  is  independent  of  the  specific  values  of  Zj  and  Z\;  and 


(ii)  The  rough  surface  z  =  Z(x)  Is  a  normally  distributed  stationary 
random  process  described  by 


W(E,)w'(z.:)  =  -L 


_ l_ 


(A. 38) 


where  ^  and  T  are  the  rms  surface  height  and  slope.  With  these  two  assump¬ 
tions,  Wagner  derived  from  Eq.  (A. 37)  an  explicit  expression  for  the  shadowing 
function  as 


H(0.)  * 

4B 


(A.39) 


where 


A 

3 


_  A.4;iA 

a 


(A. 40) 


In  the  limiting  case  A  «  1,  or  low  grazing  angles,  it  was  shown  that 


H(9,)=  JxA 


(A.41) 


For  a  general  scattering  geometry,  the  shadowing  function  in  Eq.  (A.39) 
can  be  generalized  to  describe  the  shadowing  effect  on  both  incident  and 
scattered  rays.  The  results  shown  by  Wagner  are  dependent  on  the  scattering 
angle  dg  as  follows: 


(i)  0  ^  e,  ^  r/z : 

-  4“^.  +  h) -  ^ 
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where  the  A's  and  B's  are  defined  as  in  Eq.  (A. 40); 


(ii) 


(iii) 


1C/a  ^  ^  TC-  6^  •• 


T-  a-  ^  e<  ^  Tc 

H(6J  . 


(A. 43) 


(A. 44) 
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APPENDIX  B 


Small  Perturbation  Method 


In  Appendix  A,  we  have  reviewed  several  major  solutions  that  are  based 
on  the  Kirchhoff  approximation.  The  methodology  commonly  employed  in  these 
solutions  is  usually  referred  to  as  the  tangent  plane  method.  Another 
equally  significant  class  of  solutions  that  are  constructed  or.  the  basis  of 
the  small-waveheight  approximation  essentially  follow  the  perturbation 
approach.  According  to  this  perturbation  methodology,  the  pressure-release 
boundary  condition  at  the  randomly  rough  surface  is  transformed  to  a  set  of 
boundary  conditions  at  the  unperturbed  surface  as  the  scattered  sound  field 
is  expressed  as  a  sum  of  perturbation  fields  of  various  orders.  Under  these 
newly  transformed  boundary  conditions,  the  scattered  field  is  then  solved 
either  by  the  Fourier  transform  technique  or  through  the  Helmholtz  Integral. 

This  approach  was  originated  by  Rayleigh  [CllJ  who  employed  the  small- 
waveheight  approximation  to  obtain  a  scattering  solution  for  a  periodically 
rough  surface.  Indeed,  Rayleigh's  solution  did  not  use  the  perturbation  pro¬ 
cedure  as  such  but  was  associated  with  this  methodology  because  of  its  simi¬ 
lar  approximate  boundary  condition.  Later,  Marsh  LC13-C16J  attempted  to 
obtain  an  exact  solution  for  the  problem  of  scattering  from  a  frozen  randomly 
rough  surface  through  the  use  of  Wiener's  generalized  harmonic  analysis.  Howr 
ever,  the  useful,  explicit  form  of  his  solution  was  an  approximate  one,  which 
turned  out  to  be  limited  by  the  small-waveheight  approximation.  The  rigorous 
perturbation  procedure  was  eventually  applied  to  the  surface  scattering  prob¬ 
lem  by  a  group  of  authors  (Brekhovskikh  lA3J,  Bass  and  Fuks  1A4J ,  Kuperman 
IC6-C8J ,  and  Harper  and  Labianca  [C1-C5J).  Of  these  authors,  on  the  basis  of 
their  solution  procedures,  they  were  seen  to  follow  two  distinct  approaches 
and  thus  could  be  separated  into  two  groups.  The  first  group,  including 
Brekhovskikh,  Bass  and  Fuks,  and  Kuperman,  described  the  sound  field  as  con¬ 
sisting  of  a  mean  field  and  a  perturbation  field,  and  transformed  the  pressure- 
release  boundary  condition  through  a  Taylor  series  expansion  to  a  set  of 
boundary  conditions  at  the  mean  surface.  Brekhovskikh  and  Kuperman  then  went 
on  to  solve  the  problem  by  the  Fourier  transform  technique,  which  was  also 
known  as  the  spectral  representation  analysis;  while  Bass  and  Fuks  employed 
the  Green  function  technique  to  solve  the  Helmholtz  Integral.  The  second 
group,  which  was  pioneered  by  Harper  and  Labianca,  applied  the  perturbation 
procedure  to  the  total  field  up  to  the  second  order,  transformed  the  pressure- 
release  boundary  condition  in  a  similar  fashion,  and  solved  for  the  scattered 
field  within  the  framework  of  a  normal-mode  representation.  This  last  step 
of  Harper  and  Labianca 's  approach  implies  that  they  have  recast  the  scatter¬ 
ing  problem  into  a  propagation  one,  although  not  as  completely  as  later 
authors  (J1-J6J  who  obviously  took  the  lead  from  them.  In  this  section,  we 
will  therefore  focus  on  the  works  of  Marsh,  Brekhovskikh,  and  Kuperman. 
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Bl.  Marsh 


Consider  an  incident  plane  wave  and  the  scattered  field  pg.  The 
scattered  field  is  then  expressed  by  a  plane  wave  representation  of  the  form 


r.(^)  - 


(B.l) 


where  G(ksx»ksy)  is  the  generalized  spectrum  of  the  scattered  field.  The 
pressure-release  boundary  condition  at  the  rough  surface,  Pi  +  Ps  “  0»  can 
be  written  as 


r  - 

-e 


-ik.f 


(B.2) 


Next,  the  generalized  spectrum  G  is  represented  as  a  power  series  in  (T  , 
the  rms  surface  height. 


(B.3) 


Substitution  of  Eq.  (B.3)  in  to  Eq.  (B.2)  and  collection  of  terms  of  the 
same  power  in  0“  yield  an  infinite  set  of  linear  differential-integral  equa¬ 
tions  that  in  principle  can  be  solved  for  the  coefficients  A„'s.  Each  of 
these  linear  equations  is  of  the  form 

-o* 


in  which  Z  is  the  normalized  surface  waveheight 


Harsh  IC15J  solved  a  set  of  equations  represented  by  Eq.  (B.4)  for 
m  *  0,  1,  and  2  for  dA(},  dAx ,  and  dA2.  Through  Eqs.  (B.l)  and  (B.3),  the 
covariance  function  ^  of  the  scattered  field  is  established  as 

where  is  the  spatial  covariance  function  and  Sc  the  spatial  spectrum  of 
the  surface  waves.  The  Fourier  transform  of  then  provides  the  scatter¬ 

ing  kernel 


(B.5) 


(B.6) 
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where  K  =  k.g  -  as  before  and  R  representing  the  specular  reflection  inten¬ 
sity  is  given  by 

R  =  l_  S  «r,.  (4<„-  •  (B.7) 

In  Eq.  (B.6),  the  scattering  kernel  M(ksk»  •'■sy^  describes  the  scattered 
intensity  in  the  direction  of  kg. 


B2.  Brekhovskikh 

As  indicated  earlier,  there  are  several  authors  who  essentially  followed 
a  similar  perturbation  approach.  To  save  space  and  for  the  ease  of  compara¬ 
tive  analysis,  it  is  desirable  to  present  their  results  in  a  common  terminol¬ 
ogy  despite  the  fact  that  each  author  focused  on  different  aspects  of  the 
problem.  Differences  or  similarities  will  be  pointed  out  wherever  appropri¬ 
ate.  Brekhovskikh 's  analysis  tA3,  Chapter  9j  of  the  problem  of  scattering 
from  the  ocean  surface  is  probably  the  most  succinct  presentation  and  will 
therefore  be  discussed  first. 

B21.  First-order  Perturbation  Field 

The  pressure-release  boundary  condition  for  the  total  field  p(r)  at  the 
rough  surface  z  >  Z(x,y)  is  expanded  into  a  Taylor  series  up  to  the  first 
order  in  Z, 

+  =  0  «t:  KsO  .  (B.8) 

Next,  the  total  sound  field  in  the  half  space  z  >  0  is  written  according 
to  the  perturbation  procedure. 


^  (B.9) 

where  is  the  sound  field  for  the  flat  mean  surface  and  u  is  the  first-order 
perturbation  of  the  scattered  field. 


Substituting  Eq.  (B.9)  in  Eq.  (B.8)  and  equating  terms  of  the  same  order 
of  Z*  and  Z,  the  boundary  condition  at  z  ■  Z(x,y)  becomes  the  boundary  condi¬ 
tions  at  z  *  0, 


mCk,  j) 


0. 


sue?) 


«-0 


(B.IO) 
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by 


For  an  incident  plane  wave,  the  solution  for  a  flat  interface  is  given 


( 


(B.ll) 


From  Eq.  (B.IO),  the  perturbation  u  at  z  =  0  is  then 


(B.12) 


Describing  the  surface  waveheight  Z(rj^ 
integral. 


Z(rJ  =  ^  JlK  ^ 


in  Eq.  (B.12)  by  a  Fourier 


(B.13) 


where  1C  is  the  surface  wave  vector  and  A(1C)  is  the  complex  amplitude  of  the 
constituent  Fourier  components,  one  can  obtain  a  spectral  representation  for 
the  perturbation  field  in  the  half-space  z  >  0, 


mo 

\  A(Jc)  e 


ie 

X  .  T* 


JlK 


(B.14) 


a-t  'A 


where  K.}  J  • 


Equation  (B.14)  Implies  that  within  the  framework  of  the  first-order  per¬ 
turbation  theory,  each  Fourier  component  of  the  surface  roughness  gives  rise 
to  its  own  scattered  plane  wave  in  the  direction  kg  according  to  Bragg's  law 


(B.15) 


B22.  Scattering  Cross-section  in  the  Far  Field 

From  Eq.  (B.14),  the  scattered  sound  intensity  Ig  *  <u  u*>  is  written  as 

-•0 

where  ■*  •  spectral  amplitudes  in  the  above  equation 

are  delta-correlated  [A3j,  or 


<A(ft)A*(i?)>  =  S.(tC)  JCk-Ic')  , 

S  / 
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in  which,  as  before,  SgClO  is  the  spatial  spectrum  of  the  surface  roughness. 


Substitution  of  Eq.  (B.17)  into  Eq.  (B.16)  yields 


Is  -  K  i  . 


(B.18) 


In  the  far  field  zone,  the  scattered  intensity  can  also  be  formally 
described  by 


Is  -  i'^ss(®s,'?s)'-s‘«i^,„  ;  (».19) 

where  «■«(».  A)  is  the  scattering  cross-section  per  unit  solid  angle  per  unit 
area  for  unit  Incident  intensity,  which  is  described  here  as  a  function  of 
incident  angle  0g  and  azimuthal  angle  <Pg  of  the  scattering  direction  and  is 
Integrated  over  the  flat  Insonlfled  area  »  It  was  demonstrated  by 

Br_ekhov^kikh  [A3J  that  #1^  « J ‘^5  in  the  far-field  approximation  and 
«lic  •  ***.  These  results  are  inserted  into  Eqs.  (B.18)  and 

(B.19),  which  are  then  equated  to  result  in 

nri/a. 


(B.20) 


A-  .A 


(B.21) 


in  which  kg  and  kj^  are  related  by 


(B.22) 


In  summary,  Eq.  (B.21)  describes  the  coherent  scattering  cross-section 
for  a  static  randomly  rough  ocean  surface,  which  is  seen  from  Eq.  (B.22)  to 
follow  the  familiar  Bragg's  law  of  scattering  up  to  the  first  order  of 
perturbation. 


B23.  Far  Field  Scattering  Cross-section  for  a  Moving  Ocean  Surface 

So  far  the  temporal  variation  of  the  ocean  surface  has  not  been  Included 
in  the  analysis  of  the  scattered  field.  For  this  purpose,  the  moving  surface 
displacement  should  be  represented  as  a  superposition  of  progressive  surface 

^(1.0=  (A^(K,)<e  oLK? 

'  (B.23) 


where  the  frequency  of  the  surface  wave  and  Aq(K)  is  its  real 

amplitude.  This  amplitude  is  related  to  the  spatial  spectrum  of  the  surface 
waves  So(IC)  by  the  following  relationships 

<a.(±k)a;;(±tc')>  =  ,,, 

<A.(±K)A*.(i:K')>=  0. 

^  Note  that  SqCK)  is  not  the  spatial  spectrum  of  the  motionless  sea  surface 
Ss(IO  and  is  related  to  the  latter  by 


(B.25) 


As  for  Eq.  (B.14),  the  perturbation  scattered  field  u(rg,t)  can  be 
obtained  in  a  similar  manner,  adding  the  time-dependent  factors, 

**  •i«o 

where,  with  c  being  the  sound  velocity  in  water, 

% 


JHc 


(B.26) 


(B.27a) 


Applying  the  narrow-band  approximation,  jQ,«  CO^,  one  rewrites  4^^  as 

(B.27b) 


As  K  is  replaced  by  -1C  in  the  second  integral  of  Eq.  (B.26>,  the  scat- 

] 


tered  perturbation  field  now  becomes 

I"*  x\  'I  T  Pa  ”^(‘**0+-^)'^  A  /  I?\ 


(B.28) 


where  k  *  k  •  From  Eq.  (B,28),_^the  temporal  correlation  function  of  the 
scatterl  i  fiefl,  ,  is  obtained  through  the  use  of 

relationships  given  in  Eq.  (B.24), 

*1^  -«A 


In  consideration  of  the  far  field  zone  and  a  finite  scattering  surface, 
Eq.  (B.29)  leads  to 

(r(iS^a))=  ik  k  ^ 
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(B.30) 


following  the  steps  similar  to  those  taken  from  Eq.  (B.18)  to  Eq.  (B.21). 
Equation  (B.30)  describes  two  spectral  lines  shifted  relative  to  by  the 
value  i:  Through  the  relationship  presented  in  Eq.  (B.25),  the  result¬ 

ing  scattering  cross-section  for  a  moving  surface  is  reducible  to  that  for  a 
frozen  surface  after  an  integration  of  Eq.  (B.30)  over  60  is  performed. 


B3.  Kuperman 

From  the  fundamental  concepts  initiated  by  Brekhovskikh  lA3]  and  Bass 
and  Fuks  lA4J,  Kuperman  carried  out  a  more  systematic  analysis  of  the  scatter¬ 
ing  problem  from  a  general  two-fluid  interface  using  the  perturbation  method. 
However,  his  detailed  treatment  was  mainly  focused  on  scattering  in  the  specu¬ 
lar  direction.  In  the  following,  we  will  summarize  the  major  results  contrib¬ 
uted  by  Kuperman  regarding  scattering  from  the  ocean  surface. 


B31.  Generalized  Boundary  Conditions  Transformed  by  the  Perturbation  Method 

The  geometry  of  a  general  two-fluid  interface  was  considered  with  sub¬ 
script  1  referring  to  the  incident  half-space  and  subscript  2  to  the  other 
half-space.  LetTJ^  be  the  velocity  potentials  in  the  two  half-spaces,  the 
conditions  that  velocity  and  pressure  be  continuous  across  the  interface 
could  be  stated  as 


(B.31) 


at  z  ■  Z(^).  These  conditions  were  expanded  in  a  Taylor  series  about  z  =  0 
up  to  terms  of  the  orders  and  lC8j, 


(B.32) 


m 

"iSL  I2.O 

where  one  has  used 


IjjsO  ^  *i=0 


(B.33) 


(B.34) 
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The  velocity  potential  was  decomposed  into  the  mean  field  <u.>  and 
the  stochastic  field  u^, 


"U  =  <  ■cr,>  -I-  u  j  i=  I, z  . 


(B.35) 


Substituting  Eq.  (8.35)  into  Eqs.  (B.32)  and  (B.33),  carrying  out  the 
averaging  operation,  and  a  simple  algebraic  manipulation  lead  to  two  separate 
sets  of  transformed  boundary  conditions  on  the  plane  z  =  0  for  the  mean  and 
perturbation  fields: 

(i)  For  ^U!^up  to  the  T?-  order. 


(ii)  For  u^  up  to  the  Z  order, 

f.'ii-p."!  =  <3-,  <‘oi>)  , 

Is  =  • 


(B.36a) 


(B.36b) 


(B.37a) 


(B.37b) 


Detailed  expressions  for  Fj,  F2,  Gj,  and  G2  are  given  In  Reference  lC8J. 
In  general,  the  method  of  solution  suggested  by  Kuperman  was  to  solve  for  the 
perturbation  fields  from  Eq.  (B.37)  by  the  Fourier  transform  technique  or 
spectral  analysis  and  then  to  obtain  the  mean  fields  at  the  surface  from 
Eq.  (B.36).  A  critical  step  employed  by  Kuperman  In  obtaining  a  formal  solu¬ 
tion  for  the  mean  fields  was  to  replace  the  mean  fields  In  Eq.  (B.36)  by 
solutions  to  the  flat-interface  scattering  problem.  On  one  hand,  this  greatly 
reduced  the  algebraic  complexity  but,  on  the  other  hand,  limited  the  accuracy 
of  the  resulting  mean  fields  only  to  a  lower  order  Z,  Instead  of  T?-  as  Implied 
In  the  derivation  procedure  leading  to  Eq.  (B.36). 


B32.  Reflection  Coefficient  for  Scattering  from  the  Ocean  Surface 

For  the  pressure-release  case,  the  boundary  conditions  represented  by 
Eqs.  (B.36a)  and  (B.37a)  reduce  to 


'J2=0 


0  ^ 


(B.38) 


(B.39) 
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in  which  the  unsubscripted  fields  refer  to  scattered  fields  in  medium  1.  Fol¬ 
lowing  the  procedure  briefly  described  above,  Kuperman  derived  the  expression 
for  the  mean  field  at  the  surface 


<V>^H  “  \  4  5^(k)  JIK 


(B.40) 


which  has  been  translated  into  our  present  terminology  from  Eq.  (61)  of  Ref¬ 
erence  (C8J .  Note  that  the  spatial  spectrum  of  the  surface  roughness  defined 
by  Kuperman  differs  from  Sg(K)  by  a  factor  2TC.  When  the  mean  field  is  cast 
in  the  form  of  a  flat-interface  solution 


TJ(rj)  -  e  -‘-(c.  ^  -h  R.^  **)  ^ 


(B.41) 


with  R  being  the  specular  reflection  coefficient,  then  R  can  be  identified 
from  Eq.  (B.40)  as 


R=-I+J.i.  ^4  SiK)JLK 

1  SJt  s  / 


(B.42) 


where  kg^  has  been  determined  from  Eq.  (B,27b). 

At  this  point,  it  Is  interesting  point  out  that  Brekhovskikh  1A3,  Sec¬ 
tion  9.6J  arrived  at  the  same  result  by  employing  a  rather  heuristic  approach. 
On  the  basis  of  energy  conservation,  Brekhovskikh  reasoned  that  "the  vertical 
energy  flux  through  a  unit  surface  in  the  incident  wave  had  to  be  equal  to 
the  sum  of  the  flux  in  the  reflected  and  scattered  waves.”  From  Eq.  (B.18), 
this  was  translated  into  the  following  mathematical  statement 


V  *  •  4  4e  ^  ^  » 


(B.43) 


where  and  Sg  are  the  incident  and  scattering  angles,  respectively, 
tlon  (B.43)  was  rewritten  as 


ll?r=  I- jlH 


Equa- 


(B.44) 


in  which  cos  6g  has  been  replaced  by 


c»«  9-  -  . 

’  k. 


(B.45) 
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Under  the  approximation  that  the  scattering  contribution  represented  by 
the  integral  in  Eq.  (B.44)  is  much  smaller  than  unity,  Brekhovskikh  described 
the  reflection  coefficient  as 


)l^|  ^  S  Jllc  , 

ia  SB  *  '  ) 


(B.46) 


which  is  equivalent  to  Eq.  (B.42). 

An  approximate  expression  for  |R|  was  a^so  obtained  by  Brekhovskikh  under 
the  assumption  that  the  spatial  spectrum  Sg(IC)  is  not  significant  when  K  > 
with  A, being  the  correlation  length  of  the  rough  surface.  First,  the  radical 
factor  in  the  Integrand  in  Eq.  (B.46)  was  rewritten  as 


(B.47) 


where  (©i,  described  the  direction  of  the  incident  wave  and  1C  =  (K  cose<, 
ICsina^).  The  above  assumption  implies  that  1  ,  thus  simplify¬ 
ing  Eq.  (B.47),  *  " 


56. 


(B.48) 


Hence  the  reflection  coefficient  of  Eq.  (B.46)  became 

^  T 

I I  2:  I  -  i  ^  ^  K  iK  A  oC 

6  -T 


(B.49) 


or 


I  R|  -  JT-il:* 


*  ) 


(B.50) 


where  (T 


B33.  Scattering  Kernel  and  Reflection  Coefficient  for  a  Moving  Ocean  Surface 

In  addition  to  the  above  analyses  for  a  frozen  rough  surface,  Kuperman 
also  investigated  lC6J  the  case  of  a  dynamic  ocean  surface.  In  this  extended 
case,  he  generalized  the  expression  described  by  Eq.  (B.42)  for  the  reflection 
coefficient  and  also  derived  the  scattering  kernel  for  the  nonspecular  direc¬ 
tion.  In  deriving  the  scattering  kernel,  Kuperman  employed  a  spectral  analy¬ 
sis  approach  with  the  narrow-band  approximation  essentially  similar  to 
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Brekhovskikh' s  treatment  presented  earlier,  except  that  he  focused  on  the 
scattering  kernel  but  not  on  a  far-field  scattering  cross-section  like 
Brekhovskikh.  Because  of  this,  his  scattering  kernel  is  equivalent,  except 
for  a  z-dependent  factor,  to  that  implied  by  Eq.  (B.30),  which  was  derived  by 
Brekhovskikh  for  the  far-field  case.  This  scattering  kernel  in  our  present 
notation  is  of  the  form  (see  Eq.  (8a)  of  Reference  lC8J), 


S.(-K)S'Ced-H  +  -a)]  ^ 


(B.51) 


where 


II 

(  ^1 

»5 

wJl 

0 

•■i 

•■5 

t«| 

with  ks2  given 

by  Eq. 

(B.27a). 

Note  that 

in 

the 

(B.52) 


kernel  presented  in  Eq.  (B.51)  is  Identical  to  the  scattering  kernel  derivable 
from  Eq.  (B.30). 


Inclusion  of  the  temporal  variation  also  enabled  Kuperman  to  modify  his 
expression  (Eq.  (B.42))  for  the  reflection  coefficient  to  decompose  the 
reflected  field  in  two  Doppler-shifted  components.  Starting  from  the  boundary 
conditions  represented  by  Eqs.  (B.38)  and  (£.39),  using  the  spectral  analysis 
approach,  and  replacing  the  mean  fields  with  the  flat-surface  solutions, 
Kuperman  derived  the  following  reflection  coefficient  for  a  dynamic  ocean 
surface  lC6,  Eq.  (9c) J 


(B.53) 
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APPENDIX  C 


Composite-Roughness  Approach 

In  Appendices  A  and  B,  we  have  discussed  two  fundamental  approaches  to 
the  theoretical  solution  of  the  problem  of  scattering  from  the  ocean  surface: 
the  perturbation  method  and  the  tangent  plane  method.  By  their  basic  assump¬ 
tions,  the  perturbation  method  tends  to  ignore  the  influence  of  the  large- 
scale  roughness  while  the  tangent  plane  method  could  not  reflect  the  contri¬ 
bution  of  the  small-scale  surface  waves.  A  modification  to  the  perturbation 
method  was  suggested  to  include  scattering  from  the  entire  surface  spectrum 
and  was  normally  referred  to  as  the  composite-roughness  approach.  In  the 
perturbation  method,  the  randomly  rough  surface  was  conceived  as  a  perturba¬ 
tion  from  the  flat  Interface  between  the  two  media.  By  this  principle,  the 
perturbation  method  solution  was  constructed  from  two  fundamental  steps: 
first,  the  pressure  release  condition  at  the  ocean  surface  was  expanded  about 
the  flat  interface  with  only  terms  of  low  orders  being  kept  for  consideration; 
and  second,  the  scattered  field  was  described  by  two  components:  the  pertur¬ 
bation  field  and  the  unperturbed  field,  of  which  the  latter  was  a  solution  to 
the  flat-surface  problem.  In  order  to  improve  the  perturbation  method  solu¬ 
tion,  the  composite-roughness  approach  divides  the  ocean  surface  spectrum 
into  two  parts,  high-frequency  and  low-frequency  sections.  The  low-frequency 
part  describes  the  spectral  power  of  the  large-scale  roughness  and  the  high-  • 
frequency  part  corresponds  to  the  small-scale  roughness.  The  perturbation 
procedure  is  then  applied  in  a  similar  manner  except  that  the  flat  Interface 
as  the  unperturbed  boundary  is  now  replaced  with  the  large-scale  randomly 
rough  surface.  Several  authors  have  employed  this  approach  to  investigate 
the  surface  scattering  problem  lDl-03,A4J.  In  this  chapter,  we  will  summarize 
a  heuristic  formulation  by  Bachmann  [Dl]  and  a  more  rigorous  treatment  by 
McDaniel  and  Gorman  lD3J  for  the  backscatterlng  case.  We  will  also  extend 
the  results  of  McDaniel  and  Gorman  to  the  general  scattering  case. 


Cl.  Bachmann 

Bachmann  reasoned  that  the  first-order  perturbation  theory  was  valid 
only  for  the  small-scale  rough  surface  ^5.  The  large-scale  rough  surface^L 
was  then  taken  Into  account  by  modulation  of  the  local  grazing  angle.  On  the 
basis  of  this  heuristic  conception,  Bachmann  set  out  to  Improve  the  result  of 
the  perturbation  theory  as  follows.  In  fact,  the  following  discussion  Is 
slightly  more  general  than  that  presented  In  Reference  IDIJ  by  Bachman  but 
essentially  represents  the  same  results.  Figure  Cl  Illustrates  that,  due  to 
the  presence  of  the  large-scale  roughness,  the  nominal  grazing  angle  0i  in 
the  solution  of  the  perturbation  method  should  be  replaced  by  the  local  graz¬ 
ing  •'gle  Sj  ,  which  deviates  from  61  by  the  local  slope  S  of  the  large-scale 
surface.  Therefore,  In  the  backscatterlng  case,  the  scattering  cross-section 
should  be  written  as  (see  Eq.(B.21),  for  example) 

slii.cse.  ,0)  , 

Sc  I  *■  JS'4,  *■ )  '  ) 
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where  S  is  the  two-dimensional  spatial  surface  spectrum.  The  new  backscat- 
tering  cross-section  is  then  determined  by  averaging  6^)  with  respect  to 
the  local  slope  6  , 

4,  ^*4  - 

Noting  that  ft  &  as  d.-vO  ,  Bachmann  used  an  approximation  that 
allowed  the  averaging  process  to  be  operated  over  0^  instead  of  over 
As  a  result,  this  is  formally  equivalent  to  attributing  an  inclination  to  the 
mean  sea  surface, 

+  ,  Cc.3) 

where  A&  is  determined  from  the  averaging  operation, 

4L  a  2.  ^  (C*^) 

Under  the  assumption  that  the  local  surface  slope  is  normally  distributed 
and  its  mean  value  is  zero,  the  expectation  values  of  all  odd  moments  vanish 
and  those  of  even  moments  are  given  by 

where  ^rms  i®  value  of€-.  Thus,  according  to  this  heuristic 

composite-roughness  approach,  the  backscattering  cross-section  is  described 
by 

,  O)  j  (C.6) 

in  which  6B  is  specified  by 

A6  =  5 .  (C.7) 

C.2  McDaniel  and  Gorman 

In  the  preceding  section  was  presented  an  ad-hoc  modification  to  the 
backs cat te ring  solution  of  the  perturbation  theory  employing  the  composite¬ 
roughness  concept.  On  the  basis  of  the  same  physical  picture,  McDaniel  and 
Gorman  lD3],  however,  rigorously  applied  a  perturbation  procedure  to  the 
backseat tering  case.  They  followed  the  perturbation  procedure  initiated  by 
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Labianca  and  Harper  IC1-C5J  but  considered  the  large-scale  rough  surface  as 
the  unperturbed  surface  Instead  of  the  flat  Interface.  McDaniel  and  Gorman 
started  from  a  Green's  function  form  of  the  solution  to  the  Helmhotz  equation 


(C.8) 


where  is  the^total  scattered  field,  is  the  Green  function 

to  be  chosen  appropriately  with  the  boundary  conditions,  n  Is  the  unit  vector 
normal  to  the  large-scale  surface,  and  the  Integration  Is  to  be  carried  over 
the  large-scale  surface  A-_.  Tne  scattered  field  Is  next  expressed  as  a  series 
In  the  small  parameter  ^  ,  the  rms  wavehelght  of  the  small-scale  surface^  , 
up  to  the  second  order  ® 


since  the  Integration  In  Eq.(C.8)  Is  over^l^,  the  pressure-release 
boundary  condition  at  needs  to  be  transformed  to  ^  .  This  can  be 
achieved  by  first  expanding  about  as  ^ 


(C.9) 


[U]  =  , 


(C.IO) 


where  Is  the  wavehelght  of  with  respect  to  .  Applying  the  bound¬ 

ary  condition  r\/3  wO  «  using  Eq.(C«9),  and  equating  terms  of  the  same  powers 
of  T yield  ^ 


[U.l.  =  0  , 


[U.]  =  ,  -i 

U 

[UJ  .  .  [.J.7U.]  -  [  . 


(C.lla) 


(C.llb) 


(C.llc) 


Substituting  Eq8.(C.9)  and  (C.ll)  into  Eq.(C.8)  and  equating  terms  of 
the  same  powers  of  (T  lead  to 


[&.(??) 5. 7tr (?,?.)]  j£ 

rTJ, (?.,?.)  = 

\  t  tT  *  *^|A  t-iV*!,,/.*  ♦O  i 


(C.12a) 


(C.12b) 


(C.12c) 
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The  above  set  of  integral  equations  define  the  formal  solution  of  the 
scattered  field  up  to  the  second  perturbation  order.  Since  results  obtained 
by  the  perturbation  method  indicated  that  information  on  the  scattered  field 
was  contained  in  the  first  and  higher-order  perturbation  components,  McDaniel 
and  Gorman  proceeded  to  solve  for  u.  by  talcing  two  key  steps:  (i)  The  Green's 
function  was  chosen  to  be  the  Green's  function  for  the  problem  of  scattering 
from  the  large-scale  surface  alone,  that  is  Gq  vanishes  on  and  (ii)  The 

Kirchhoff's  approximation  was  assumed  for  the  unperturbed  component.  From 
Eqs.(C.ila)  and  (G.12a),  the  first  step  implies 


In  addition,  the  Kirchhoff's  assumption  on  V.  leads  to 


(C.13) 


(C.14) 


where  V.  is  the  incident  acoustic  field. 

me. 

For  a  point  source  and  a  receiver  located  in  the  far  field,  one  can 


write 


I'll 


and  obtain 


(C.15) 


=  aii.r  _ 

Substitution  of  Eqs.(C.13)  and  (C.16)  into  Eq.(C.12b)  yields 

.mU) t ..  Wr,;.t.-,-(CC  vi} 


(C.16) 


•  (C.17) 


From  the  scattered  field  described  by  the  above  equation,  McDaniel  and 
Gorman  performed  the  ensemble  averaging  operation  over  the  large-scale  sur¬ 
face  to  obtain  the  scattered  Intensity 


=  -^<(“  OVXfj>  , 


(C.18) 


where  A  is  the  insonified  area  and 


K  -  . 


(C.19) 
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From  Eq.  (C.18)  the  general  scattering  cross-section  is  given  by 


which  is  reduced  to  the  backscattering  case  as 


(C.20) 


<r^.  4<C^-C)‘'^  SCAic«s6.jO)  ^ 

which  is  identical  to  the  earlier  result  (compare  Eq.  (C.l))  obtained  by  the 
heuristic  formulation. 


C.3  Two-dimensional  Scattering  Cross-Section 

Above  one  has  seen  the  derivation  of  the  scattering  cross-section  by  the 
composite-roughness  approach  in  both  heuristic  and  rigorous  manners.  But 
only  for  the  backscatterlng  case,  this  cross-section  could  be  reduced  to  a 
practically  useful  form  (see  Eqs.  (C.6)  and  (C.7))  through  the  use  of  an 
approximate  averaging  procedure  suggested  by  Bachmann.  In  this  section,  we 
will  employ  a  similar  procedure  to  perform  the  averaging  operation  required 
in  Eq.  (C.2U)  for  the  case  of  two-dimensional  scattering,  or  scattering  in 
the  plane  of  incidence.  In  this  case,  one  can  write  the  term  that  needs  to 
be  averaged  in  Eq.  (0.20)  as 


<(5 


(C.22) 


As  in  Section  C.l,  the  scattering  cross-section  averaged  over  the  local 
slope  can  be  formally  represented  by 


(C.23) 


in  which  the  average  inclination  Ae  is  determined  by  the  following  averaging 
operation 


(C.24) 


Again,  assuming  the  same  statistical  description  for  the  local  slope  G 
as  in  Section  C.l,  one  can  write  Eq.  (C.24)  as 


(C.25) 
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which  can  be  solved  for  A6 


)i  J 


where 


+  ;  P«=  . 

Notice  that  Eq.  (C.26)  is  reducible  to  Eq.  (C.7)  when  • 

the  two-dimensional  scattering  cross-section  is  finally  of  the  form 

with  AB  being  given  in  Eq.  (C.26). 


(C.26) 

) 

(C.27) 

Therefore, 

(C.28) 
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APPENDIX  D 


Empirical  Modeling 

So  far  we  have  been  concerned  with  descriptions  of  ocean  surface  scatter 
ing  which  are  based  on  cohesive  theoretical  formulations.  There  exist,  how¬ 
ever,  other  descriptions  of  the  surface  scattering  cross-section  which  essen¬ 
tially  represent  some  forms  of  empirical  fitting  of  data  obtained  in  field 
measurements.  In  view  of  the  considerably  complex  nature  of  the  surface  scat 
taring  problem  and  due  to  limited  validity  of  most  existing  theories,  these 
empirically-derived  descriptions  under  appropriate  conditions  have  proved  to 
be  rather  useful.  A  major  constraint  in  their  use  has  to  do  with  the  fact 
that  they  are  available  only  for  the  backseat tering  case.  In  this  appendix, 
we  will  discuss  two  of  the  more  established  empirical  models  of  surface  back- 
scattering  that  are  mentioned  most  often  in  the  literature. 


D.I  Chapman-Harris 


The  Chapman-Harris  [E4]  empirical  model  of  surface  backscattering  repre¬ 
sents  a  best  fit  of  measurements  made  in  the  early  sixties  for  a  range  of 
wind  speeds,  grazing  angles,  and  frequencies.  In  these  data,  backscattering 
strength  is  given  in  octave  bands  from  400  Hz  to  6.4  KHz  for  0-40°  grazing 
angle  and  0-30  knot  wind  speed.  For  most  reported  frequency  bands  and  wind 
speeds,  these  backscattering  data  when  plotted  against  grazing  angle  are  char 
acterlzed  by  two  distinct  features:  first.  In  the  range  10-40°,  backscattter 
Ing  strength  Is  strongly  dependent  on  grazing  angle;  second,  at  lower  grazing 
angles  backscattering  strength  Is  relatively  Independent  of  grazing  angle. 
Intuitively,  Chapman  and  Harris  attributed  these  two  distinct  scattering  fea¬ 
tures  to  contributions  from  two  different  mechanisms,  the  first  due  to  scat¬ 
tering  from  surface  roughness  and  the  second  due  to  volume  scattering  from  a 
layer  of  Isotropic  scatterers.  Employing  only  those  portions  of  the  scatter¬ 
ing  strength  versus  grazing  angle  curves  corresponding  to  roughness  ‘scatter¬ 
ing,  Chapman  and  Harris  constructed  an  empirical  formula  of  backscattering 
strength  as  a  function  of  grazing  angle,  acoustic  frequency,  and  wind  speed. 
This  formula  was  obtained  In  the  form 


_  4a.-4  ^  ^ 


(D.I) 


In  which 


(D.2) 


where 

=  backscattering  strength  In  dB/^Pa, 
6*  ■  grazing  angle  In  degrees, 

‘VT  >  wind  speed  In  knots,  and 
■  acoustic  frequency  In  Hz. 
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Although  the  Chapman-Harris  formula  has  been  used  extensively,  attention 
must  be  drawn  to  the  fact  that  its  domain  of  validity  should  be  constrained 
within  10-40®  grazing  angles,  0-30  knot  wind  speeds,  and  400-6400  Hz 
frequencies. 


D.2  Martin 

Martin's  approach  [ElJ  to  empirical  modeling  to  the  backscattering  prob¬ 
lem  is  slightly  different  to  that  of  Chapman  and  Harris.  Instead  of  direct 
fitting  of  experimental  measurements,  Martin  constructed  his  model  on  the 
basis  of  several  theoretical  results  and  used  available  acoustic  and  optical 
reverberation  data  to  estimate  some  key  parameters  built  in  these  theoretical 
results.  This  observation  will  become  clearer  in  the  description  of  Martin's 
model  below. 

Martin  perceived  the  surface  scattering  strength  as  being  composed  of 
three  elements;  roughness  scattering  ),  sublayer  bubble  scattering  ('T'b), 
and  facet  reflection  and  diffraction  (^)*  The  total  backscattering  strength 
is  then  written  as 


T  =  +  v;  .  (D.3) 

For  the  roughness  scattering  element,  Martin  used  the  backscattering 
cross-section  derived  by  Kuo  tC17J  as  an  extension  of  Marsh's  theory  for  the 
case  of  an  Impedance  interface.  This  result  is  indeed  reducible  to  that  of 
Marsh  or  the  first-order  perturbation  theory  for  the  case  of  a  pressure- 
release  surface.  In  our  notation,  the  roughness  scattering  component  is  then 
given  as 


(D.4) 


in  which  Martin  suggested  to  deduce  the  two-dimensional  surface  spectrum 
from  reverberation  data  if  not  otherwise  available.  Inferring  from  the  avail¬ 
able  data,  Martin  also  proposed  to  modify  the  use  of  the  above  expression 
near  normal  incidence.  Specifically,  he  suggested 

=  .  (0.5) 


With  regard  to  the  second  element  of  bubble  scattering,  Martin  applied 
the  general  theory  of  volume  scattering  developed  by  Tatarski  IV.  I.  Tatarski, 
"Wave  propagation  in  a  turbulent  medium,"  McGraw-Hill,  New  York,  81-90  (1961)J 
to  scattering  by  sublayer  bubbles.  The  general  result  of  Tatarski  was  derived 
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for  an  isotropic  distribution  of  volume  scatterers  such  as,  for  example,  tem¬ 
perature  fluctuation,  salinity  fluctuation,  or  subsurface  bubbles.  The  back- 
scattering  strength  for  these  types  of  scatterers  is  of  the  form 


6 


(D.6) 


where,  in  the  case  of  bubbles, 

^  -  characteristic  length  of  bubbles, 

effective  depth  at  which  scattering  of  k^  takes  place,  and 
=  three-dimensional  power  spectrum  of  x. 


In  absence  of  actual  data,  the  biased  power  spectrum,  which  is  the  term 
contained  in  the  curly  brackets  in  Eq.  (0.6),  was  estimated  by  Martin  from 
the  acoustic  reverberation  data.  This  results  in 
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) 


(0.7) 


in  which 

^  =  1.25  X  10“^,  a  fitting  constant, 
■=  wind  speed  in  knots,  and 
=*  incident  wave  number  in  cm”l. 


Martin  also  proposed  to  modify  the  use  of  the  above  fitting  formula, 
Eq.  (0.7),  as  dictated  by  the  involved  data  such  that 


V-  ^  S’  (^\r  e  ^-Deni’S  )  ^ 

(^V  ^  'S-A:vw^«)  =  ^ 


(0.8) 


and  to  restrict  its  application  to 


0. 0? 


^  ^  5"  c***  •-v  0.6  4  ^ 


(0.9) 


Finally,  the  contribution  by  facet  reflection/diffraction, 
own  analysis  lEZ]  ,  was  modeled  as 


based  on  his 


(0.10) 


in  which 


(0.11) 
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f  tl  fff 

where  ^  ,  O"  ,  and  O’  are  the  variances  of  the  first,  second,  and  third 
derivatives  of  the  surface  elevation.  Since  these  statistics  are  often  not 
available,  Martin  proposed  to  estimate  and  f rom  the  acoustic  and  optic 
data  of  surface  wave  spectra.  The  results  are 

X 

-h  l^)  (D.12) 

and  for  ^  at  normal  incidence 

(  ^  J- ^  “y  ^  ^0  inis 

~  I  Ab  ^  \r  >  J-a  iwis  . 

where  the  wind  speed  v  is  in  knots. 

In  summary,  one  notes  that  Martin's  model  is  a  synthesis  of  theoretical 
analyses  available  at  the  time,  which  are  then  correlated  to  acoustic  and 
optic  data  on  acoustic  reverberation  and  surface  wave  spectra  to  estimate 
the  involved  environmental  parameters.  The  final  output  is  an  analytical- 
empirical  description  of  the  backscattering  strength  as  a  function  of  grazing 
angle,  wind  speed,  and  frequency.  The  physical  rationale  to  construct  such  a 
model  is  driven  primarily  by  the  cognition  that  acoustic  scattering  from  the 
sea  surface  is  dominated  by  three  mechanisms  in  three  overlapping  ranges  of 
grazing  angles:  bubble  scattering  (O^-SO*),  roughness  scattering  (30'’-70‘’), 
and  facet  reflection  (70‘’-90®). 
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APPENDIX  E 


Empirical  Modeling  of  Surface  Wave  Spectn 

The  spectrum  of  ocean  surface  waves  is  the  most  critical  input  that 
determines  the  performance  of  any  scattering  theory.  Since  the  processes 
governing  generation  and  development  of  ocean  waves  are  very  complex,  a  cohe¬ 
sive  theory  that  is  expected  to  provide  an  adequate  description  of  these  sur¬ 
face  waves,  including  the  resultant  wave  spectra,  until  now  has  been  unavail¬ 
able.  Consequently,  most  computations  that  require  spectra  as  input  have  made 
use  of  descriptions  empirically  derived.  These  empirical  relationships  should 
be  used  with  care  due  to  lack  of  universal  agreement  and  considerable  errors 
associated  with  many  of  the  measurement  techniques. 


Of  all  the  various  types  of  wave  motion  in  the  ocean,  we  will  be  only 
concerned  with  those  surface  waves  which  are  generated  by  the  interaction  of 
the  wind  with  the  ocean  surface.  An  Important  distinction  should  be  made 
here  between  "gravity"  and  "capillary"  waves.  Surface  gravity  waves  are 
caused  by  the  wind  and  propagate  under  the  restoring  force  of  gravity.  Their 
wavelengths  range  from  about  10  cm  to  about  1  km,  with  maximum  energy  density 
typically  centered  about  the  wavelength  of  150  m.  Popular  names  for  gravi.ty 
waves  include  "sea  waves"  for  rough  irregular  waves  near  a  storm  area,  and 
"swells"  for  smooth  almost-slnusoidal  waves  at  some  distance  from  the  storm 
center.  On  the  other  hand,  capillary  waves  are  also  generated  by  the  wind 
but  propagate  under  the  restoring  force  of  surface  tension.  They  are  some¬ 
times  referred  to  as  "ripples"  with  wavelengths  less  than  about  1  cm  and 
their  amplitude  considerably  less  than  that  of  gravity  waves.  Caution  should 
be  taken  to  differentiate  between  these  two  types  of  surface  waves  because 
they  are  characterized  by  two  different  dispersion  relationships.  In  this 
study,  the  term  "ocean  surface  waves"  are  used  to  denote  both  types  and  the 
following  combined  dispersion  relationship  is  valid  for  the  entire  range  of 
wavelengths , 


A 

U)  - 


gk  (1  +  — ), 


(E.l) 


where 


is.  -  < 

in  which 

•  angular  frequency, 
k  *  wave  number, 
g  ■  gravitational  acceleration, 
P  ■  water  density,  and 
o  *  surface  tension. 
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The  first  serious  attempt  to  model  the  one-dimensional  spectrum  was 
probably  made  by  Phillips  IK.14J  for  those  frequencies  that  are  higher  than 
that  of  the  peak  in  the  observed  spectrum  and  lower  than  frequencies  for 
which  surface  tension  is  important.  This  range  of  frequencies  is  known  as 
the  "equilibrium  range",  in  which  the  stability  of  the  water  surface  is  main¬ 
tained,  i.e.,  when  the  surface  waves  have  reached  equilibrium  under  the  given 
wind  conditions  in  the  sense  that  they  cannot  grow  any  higher  without  break¬ 
ing.  On  the  basis  of  dimensional  analysis,  Phillips  derived  that  the  one¬ 
dimensional  energy  spectrum  decreased  with  increasing  frequency  (f)  as  the 
inverse  fifth  power  of  the  frequency,  or 


Sp(f )  =  ^g^f  ^  , 


(E.2) 


where >u  =  18.3  x  10~^  (MKS),  Phillips'  dimensionless  equilibrium  range  constant. 

In  the  absence  of  measurements  of  surface  wave  energy  as  a  function  of 
fetch,  wind  duration  and  water  depth,  the  modeling  efforts  continued  to  be 
developed  for  idealized  conditions.  One  of  such  considerations  indeed  led  to 
the  fairly  popular  and  useful  concept  of  a  "fully-developed  spectrum".  The 
physical  idea  is  as  follows.  In  an  ocean  initially  at  rest,  a  steady  wind 
commences  to  blow  over  an  infinite  surface.  Waves  of  high  frequency  will 
grow  fast  and  reach  equilibrium  quickly.  Waves  of  successively  lower  and 
lower  frequency  will  then  grow  to  equilibrium,  filling  Phillips'  equilibrium 
range.  The  waves  of  very  low  frequency  will  grow  to  a  certain  extent  until 
their  phase  velocity  becomes  faster  than  the  wind  speed.  One  expects  then 
both  low  and  high  cutoffs  in  frequencies  with  a  peak  at  about  the  frequency 
for  which  the  wave  phase  speed  equals  the  wind  speed.  Therefore,  after  a  suf¬ 
ficiently  long  duration,  one  anticipates  that  the  entire  spectrum  will  reach 
a  steady  state  and  will  be  fully  developed.  Although,  it  appears  limited, 
this  conceptualized  spectrum  has  proved  to  be  rather  useful  for  description 
of  surface  waves  in  the  open  ocean.  Pierson  and  Moskowitz  [K10,KllJ  proposed 
the  following  fully-developed  spectrum  model  to  fit  their  observed  data, 

e  ,  (E.3) 


where 


^  =  g/u, 

oi  =  0.0081  (MKS), 

^  =  0.74  (MKS),  and 
u  =  wind  speed  (m/s). 

This  spectrum  model  was  extensively  used  in  many  calculations  and  was 
continually  updated  with  new  measurements  for  various  frequency  ranges  lK3,K5, 
K7,K9J.  The  latest  version  developed  by  Pierson  IK4J  represents  a  systematic 
synthesis  of  several  data  sources  to  cover  the  broadest  possible  frequency 
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range.  It  is  this  version  that  will  be  used  in  the  present  study.  In  gen¬ 
eral,  Pierson  organizes  the  one-dimensional  fully-developed  spectrum  into  six 
ranges,  each  of  which  is  described  by  two  sets  of  expressions.  One  set  is  to 
express  the  spectrum  in  terms  of  frequency  and  the  other  in  terms  of  wave  num¬ 
ber.  They  are  related  to  one  another  through  a  specific  dispersion  relation¬ 
ship,  which  is  dictated  by  the  dominant  generation  mechanism,  either  gravity 
or  capillary,  for  that  particular  frequency  range.  This  model  is  summarized 
as  follows: 


A.  Pierson-Moskowitz  Range 


Frequency: 

Wave  Number:  [0,kAJ 
Dispersion  Relationship: 
Empirical  Formulae: 

9t"  - 

Slio)  =  e- 

A 


a 

=  kg  (gravity) 


a( 


Parameters; 


oC 


=  0.0081  MKS 
=  0.74  MKS 


(E.4) 

(E.5) 


B.  Stacy  Range 


Frequency:  I 
Wave  Number:  ll 

Dispersion  Relationship;  to  =  gk  (gravity) 
Empirical  Formulae: 


,  '4 


) 


(E.6) 

(E.7) 


Parameters: 


Ht/\o 

(,T( 


Tf  =  0.0543  MKS 
S  =  0.002717  MKS 
e.  =  2.53  MKS 
U«.p  =  35.8  cm/s 
=  12.0  cm/  s 

u*  =  friction  velocity  (  cm/s  ) 

minimum  friction  velocity  (  cm/s  ) 


C.  Kitaigorodskii  Kange 


Frequency: 

Wave  Number:  ^ 

Dispersion  Relationship:  w  =  kg  (gravity) 
Empirical  Formulae: 


+  ^  i 


(E.8) 


"  I  s*(«l  f '  7p 


(B.9) 


Parameters: 


ft)  : 


(410/^ 

75.76  cm/s 


D.  Leykin-Rosenberg  Range 


Frequency:  l*)^,  s!  . 

Wave  Number:  1 4(^  •v 
Dispersion  Relationship: 
Empirical  Formulae: 


*  4r 

40  =  gk  (1  +  -jx)  (gravity  and  capillary) 


%(**»)- 

F(i2)= 

1  “ 


,  +  3F*(u)/a)^) 


S(*)=  S^d,)  (!/■«'>)“’ 


(E.IO) 


(E.ll) 
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Parameters : 


•*>»  =10  71: 

363.0  !/*« 


E.  Mitsuyasu-Honda  Range 


Frequency:  IW^.WgJ 
Wave  Number: 

Dispersion  Relationship: 
Empirical  Formulae: 


A  M 

•3-  JU 

w  =  gk  (1  +-7x)  (gravity  and  capillary) 


=  (^Aft)(;LTC/a))^ 


Parameters: 

<fe '  (‘liCVJ'i 

*1  =  0.0000875  MKS 
f  =  5  -  log  u* 
i£-=  1.473  (10)"^ 


(E.12) 

(E.13) 


F .  Cox  Range 


Frequency:  t<^,»o] 

Wave  Number:  [l^ioo  J  x  )  a. 

Dispersion  Relationship:  ‘^  =  g  wK  (capillary) 
Empirical  Formulae: 


(E.14) 


(E.15) 


,  .  .  5  i  iV3 

S^(u.)  =  ^£,  “> 

SpU)  = 


It  should  be  noticed  that  in  this  model  the  only  environmental  factor 
affecting  the  resulting  spectrum  is  the  friction  velocity  u*(in  cm/s),  which 
is  directly  related  to  the  wind  speed  u  (in  m/s)  measured  at  a  certain  height 
h  (in  m)  above  the  sea  surface.  In  this  particular  case,  the  Pierson  model 
was  developed  with  the  use  of  the  well-established  Cardone  curve  IK15J,  exper¬ 
imentally  relating  u*  to  u  at  a  specified  height  h  ■  19.5  m.  This  curve  is 
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reproduced  here  in  Fig.  El.  The  following  set  of  values  of  u  and  u*  extracted 
from  the  Cardone  curve  will  be  used  in  subsequent  computations: 


u  (m/s) 

u*  (cm/s) 

3.5 

12.0 

7. A 

24.0 

12.3 

48.0 

20.6 

96.0 

In  support  of  the  modeling  effort  on  ocean  surface  scattering,  a  computer 
code  was  developed  at  NRL  to  Implement  the  above  Pierson  model  describing  the 
one-dimensional  fully-developed  spectrum  of  surface  waves  in  both  formats: 
spectral  density  (m^/Hz)  versus  frequency  (Hz)  and  spectral  density  (m^)  ver¬ 
sus  wave  number  (1/m).  Figures  E2  and  E3  illustrate  the  typical  results  for 
the  chosen  values  of  wind  speed.  Also  indicated  in  Fig.  E2  are  the  six  fre¬ 
quency  ranges  considered  by  Pierson  in  developing  this  model. 

So  far  only  wind  speed  has  been  taken  into  consideration  in  discussing 
about  the  surface  wave  spectrum.  In  fact,  other  than  the  spectral  distribu¬ 
tion  of  surface  waves,  one  would  also  like  to  learn  about  how  the  energy  is 
distributed  in  various  azimuthal  directions  given  the  knowledge  of  the  mean 
wind  direction.  This  turns  out  to  be  essential  information  if  one  desires 
to  obtain  a  more  complete  picture  of  scattering  from  the  ocean  surface.  In 
general,  the  scattering  strength  is  influenced  by  the  specific  wind  direction 
as  strongly  implied  by  all  scattering  theories.  The  primary  environmental 
input  is  always  the  two-dimensional  spectrum  of  surface  waves,  not  the  one- 
dimensional  one. 

In  reality,  much  less  is  known  about  the  two-dimensional  or  directional 
spectrum  than  the  one-dimensional  spectrum.  Present  knowledge  is  based  on  a 
very  few  attempts  to  measure  the  directional  spectrum,  resulting  in  several 
empirically-derived  descriptions  put  forward  to  summarize  these  observations. 
Two  of  the  more  comprehensive  sets  of  such  measurements  were  performed  by 
Tyler  et  al  lK8j,  who  employed  the  radio  wave  synthetic  aperture  technique, 
and  Mitsuyasu  et  al  lK6J ,  who  made  their  measurements  using  the  cloverleaf 
buoy.  In  this  effort,  we  will  use  the  azimuthal  distribution  model  proposed 
by  Mitsuyasu  et  al  in  combination  with  the  Pierson  one-dimensional  model  to 
construct  a  two-dimensional  spectrum  of  surface  waves. 

The  two-dimensional  spectrum  S2(f.6)  can  be  expressed  in  the  form 
S(4i®)  s=  S,(-f)  A(i,0)  ,  (E.16) 


where  Si(f)  is  the  one-dimensional  spectrum  and  A(f,6)  is  an  angular  distri¬ 
bution  function  with  6  denoting  the  angle  relative  with  the  wind  direction. 
The  function  A(f,6)  should  be  normalized  such  that 
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(E.17) 


AK 

^  A{-f,e)J.e  =  1  . 

0 

In  deriving  a  specific  form  for  A(f,e),  Mitsuyasu  et  al  assumed  the 
following  analytical  form  for  A(f,e) 

A(4, &)=  Nl(s)|«=«Sj|  (E.18) 

of  which  the  parameter  s  was  a  function  of  frequency  determined  by  fitting 
the  observed  data  to  this  assumed  form  and  the  function  N(s)  was  specified 
via  the  normalization  condition  of  Eq.  (E.17),  leading  to 


MCs) 


K  r(A»+i) 


where  T'  is  the  Gamma  function. 


(E.19) 


The  function  A(f,e),  as  prescribed  by  Eq.  (E.18),  represents  a  symmetric 
bell-shaped  distribution  centered  about  the  wind  direction,  the  width  of  which 
depends  solely  on  the  frequency-dependent  value  of  the  fitting  parameter  s(f). 
Figure  E4  ilustrates  this  angular  distribution  for  representative  values  of  s. 
There  it  is  seen  that  most  of  the  energy  is  strongly  concentrated  in  the  wind 
direction  for  high  values  of  a  while  it  is  more  evenly  spread  out  if  the  val¬ 
ues  of  6  are  low.  The  exact  functional  dependence  of  s  proposed  in  the 
Mitsuyasu  model  was  given  as 


i  <‘i'm 


) 


(E.20) 


where  f  is  the  dimensionless  frequency  defined  by  f  “Wu/g,  f*,  denotes  this 
quantity  corresponding  to  the  peak  frequency  of  Sj,  and  s©  “  11.5,  a  fitting 
constant.  The  above  estimate  of  the  s  parameter  is  a^reasonable  representa¬ 
tion  of  the  presented  data  except  at  the  maximum  £  ~  >  where  s  undergoes  a 

rather  sudden  change  in  value  in  the  form  of  a  very  8,har£  peak.  Here  we  pro¬ 
pose  to  remove  this  unrealistic  behavior  of  8(f)  at  £  “  £«,by  replacing  the 
fitting  expression  of  Eq.  (E.20)  with  a  slightly  different  expression  for 
8(f)  without  significantly  changing  the  original  estimate  for  s(f).  The  modi¬ 
fied  fitting  formula  is 


1  ^(S'/?r)«r [-(?/?.)');  ■ 

Both  estimates  of  8(f),  original  and  modified,  are  presented  in  Fig.  E3 
for  the  four  chosen  values  of  wind  speed.  This  demonstrates  that  the  modified 
expression  for  8(f)  is  not  substantially  different  from  the  original  one  pro¬ 
posed  by  Mitsuyasu  et  al,  making  it  realistically  well-behaved  near  the  peak 
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frequency.  Figures  E4  and  E5  combined  show  that  the  Mitsuyasu  azimuthal 
dependence  predicts  that  the  angular  distribution  will  be  relatively  well- 
focused  near  the  peak  frequency  and  become  broader  as  the  frequency  goes 
farther  away  from  it.  ' 

In  applying  the  Mitsuyasu  distribution  to  construct  a  two-dimensional 
spectrum,  we  substitute  Sj^Cf)  with  the  Pierson  one-dimensional  fully-developed 
spectrum.  This  will  result  in  a  complete  description  of  the  surface  wave 
intensity  in  both  spectral  and  angular  dimensions.  With  this  representation 
for  the  ocean  surface  waves,  we  have  the  necessary  Information  required  for 
input  into  all  surface  scattering  models  of  Interest.  A  graphical  illustra¬ 
tion  of  the  two-dimensional  spectrum  is  presented  in  Fig.  E6  for  a  typical 
value  of  wind  speed.  Therefore,  on  the  basis  of  Mitsuyasu 's  model  of  azi¬ 
muthal  dependence,  variation  of  the  surface  wave  spectrum  as  a  function  of 
the  azimuth  is  only  significant  for  a  limited  frequency  range  about  the  peak 
frequency.  The  implications  of  this  on  acoustic  surface  scattering  is  that 
the  directionality  will  be  frequency  dependent  for  a  developed  sea  and  poten¬ 
tially  vary  at  a  fixed  frequency  as  the  sea  develops.  This  shows  that 
extremely  careful  experimental  measures  will  be  required  to  test  these  com¬ 
bined  acoustic  and  surface  wave  hypotheses. 
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Fig.  2  —  Variation  of  the  isotropic  backscattering  strength  with  wind  speed  computed  by  the  first-order  pertur¬ 
bation  method  at  acoustic  frequencies  of  a)  200  Hz  and  b)  SO  Hz. 
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GRAZING  ANGLE  (Deg) 

Fig.  2  —  (Continued)  Variation  of  the  isotropic  backscttering  strength  with  wind  speed  computed  by  the  first 
order  perturbation  method  at  acoustic  frequencies  of  a)  200  Hz  and  b)  50  Hz. 


GRAZING  ANGLE  (Deg) 

Fig.  3  —  Variation  of  the  isotropic  backscattering  strength  with  wind  speed  computed  by  the  composite¬ 
roughness  approach  at  acoustic  frequency  of  200  Hz.  Different  ocean  surface  models  used  as  input  to  this  solu¬ 
tion  are  compared  between  the  Phillips  spectrum  and  the  six-range  Pierson  spectrum  adapted  in  the  present 


Fig.  S  —  The  up-shifted  and  down-shifted  components  of  the  backscattering  strength  computed  by  the 
composite-roughness  model  are  evaluated  for  four  different  wind  directions  (with  respect  to  the  incident  direc¬ 
tion  projected  on  the  surface):  a)  100  deg  (almost  cross-wind);  b)  120  deg;  c)  ISO  deg;  d)  180  deg  (up-wind). 
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Fig.  5  —  (Continued)  The  up-shifted  and  down-shifted  components  of  the  backscattering  strength  computed  by 
the  composite-roughness  model  are  evaluated  for  four  different  wind  directions  (with  respect  to  the  incident 
direction  projected  on  the  surface):  a)  100  deg  (almost  cross-wind);  b)  120  deg;  c)  150  deg;  d)  180  deg  (up¬ 
wind). 
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Fig.  5  —  (Continued)  The  up-shifted  and  down-shifted  components  of  the  backscattering  strength  computed  by 
the  composite-roughness  model  are  evaluated  for  four  different  wind  directions  (with  respect  to  the  incident 
direction  projected  on  the  surface):  a)  100  deg  (almost  cross-wind);  b)  120  deg;  c)  ISO  deg;  d)  180  deg  (up¬ 
wind). 
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Fig.  6  —  The  incoherent  (up-shifted  plus  down-shifte)  scattering  strength  as  a  function  of  the  scattering  angle 
with  the  incident  grazing  angle  being  30  deg.  Computations  are  performed  using  both  the  perturbation  and 
composite-roughness  methods  for  2(X)  Hz  frequency  and  3.5  m/s  wind  speed  under  two  different  wind  condi¬ 
tions:  cross-wind  and  up-wind. 
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Fig.  A1  —  Scattering  geometry  and  coordinate  system  employed 
in  discussions  on  the  tangent  plane  method. 
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WIND  SPEED  (knots) 


Fig.  El  —  The  Cardonc  curve  relating  the  surface  friction  velocity  and  wind  speed 
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Fig.  E2  —  Six-range  Pierson’s  model  of  one-dimensional  fully-developed 
specmun  of  ocean  surface  waves  in  the  frequency  domain. 
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Fig.  E3  —  The  same  Pierson’s  model  as  in  Fig.  E2  expressed  in  the  wave  number  domain 
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Fig.  E5  —  Empirical  estimates  of  the  s  parameters  for  a  given  wind  speed  of  7.4  m/s.  Both  Mitsuyasu’s  origi 
nal  estimate  and  our  modification  are  presented  for  comparison. 
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—  Two  dimensional  spectrum  of  surface  wave  constructed 
I’s  models  is  illustrated  for  a  typical  wind  speed  of  12.3  m/s. 
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